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ABSTRACT 

Title  of  Thesis s  Extensions  of  the  Laplace  Cascade  Method. 
John  Hilary  Billings,  Doctor  of  Philosophy,  i960. 
Thesis  directed  "by:  Research  Professor  Joaquin  B.  Diaz. 

The  Laplace  cascade  aethod  is  concerned  with  second  order  linear 
hyperbolic  equations  of  the  form 

(1)  u^  +  a(x,y)  u^  +■  h(x,y)  uy  4-  c(x,y)  u  -  0  . 

The  substitutions  u-,  ■  u_  4-  au  and  u  -  =  ux  +■  hu  lead  to  the 
equations 

(2) 


u,  +  hu^  -  hu  • 


and 


(3)  ux  +  aU;L  -  ku  -  0  , 

where  h  m   a  +  ah  -  c  and  k  »  h  +  ah  -  c  are  the  two  Darboux 

x  y 

invar iants.  If  either  h  or  k  vanishes,  (l)  has  "been  reduced  to  a 
system  of  two  first  order  equations,  while  if  neither  vanishes  the 
equation  may  he  cascaded  in  two  directions.  Solving  (2)  (or  (3))  for 
u  in  terms  of  u,(or  u_^),  and  substituting  the  resulting  expression 
into  (l)  yields  an  equation  for  UtCu^)  of  the  same  form  as  (l)  hut 
with  new  coefficients,  in  general.  This  process  may  he  iterated,  forming, 
a  chain  of  equations,  until  either  the  original  equation  reappears,  or  one 
of  the  corresponding  invariants  vanishes. 

An  extension  of  Volterra's  product  integral  to  the  non-homogeneous 
system 


n 


(*0    dx       Jsl 


u^b)  =  uiQ     |   i  =  l,2,...,n 


is  made  first.  Then  the  Laplace  Method  is  extended  to  systems  of 
second  order  hyperbolic  equations,  of  the  form 


=  0  , 


(5)  -^ +  Z! a^ s  uj + y  *«  ^ u* +  z  vj 

^xdy      J  =  l       ^  X     J*l     Jy^   J=l 

i  —  1,2,  ...,n  # 
Matrix  notation  is  promptly  introduced,  and  equation     (5)     is  rewritten 

as 

(6)  U       +    AU      +    BU+    CU     -    0 

The  two  related  matrix  invariants  for  (6)  are 

H  «  A^  +  BA  -  C 

K  =  B  4-  AB  -  C  , 

and  the  chain  of  equations  is  developed  as  for  the  single  equation.  If 
either  H  or  K    is  identically  zero,  the  resulting  two  systems  of  first 
order  equations .may  he  solved  by  employing  the  ah ove -mentioned  extension  of 
Volterra'e  product  integral. 

The  invariance  of  H  and  K  are  discussed  in  a  manner  analogous 
to  the  invariance  of  the  functions  h  and  k  of  the  single  equation. 
This  is  followed  by  a  consideration  of  periodic  systems,  i.e.  systems  such 


that  after  „j  Iterations  the  original  equation  reappeaapa.  This  dis- 
cussion results  in  two  theorems,  the  first  of  which  is 

Theorem  I  -  A  system  of  equations  of  the  form  (6)  having 
constant  matrix  coefficients  A  and  B,  can  he  reduced  to  the  form 

u   =  HB 

xy 

by  a  change  of  variables  U  =  AU  ,  if  an  only  if  AB  *  BA.  The  second 
theorem,  arising  from  systems  of  period  two  leads  to  a  discussion  of  the 
form  of  the  solution  to  the  matrix  analog  of  Liouville's  equation 


Discussion  of  the  form  of  the  solution  when  the  chain  terminates 
after  a  finite  number  of  iterations  leads  to  two  further  theorems  which 
are  completely  analogous  to  theorems  proved  hy  Darboux  for  the  single 
equation. 

The  second  extension  of  the  Laplace  cascade  method  is  made  to  the 
third  order  linear  hyperbolic  equation  in  three  independent  variables  of 
the  form 

(7)   ^ra^r  ^  +  cu^  +  di^  +  euy  +  fus  +  ga  =  0  . 

Here  the  number  of  invariant  functions  to  be  consider  jumps  from  two  to 
eighteen.  The  nature  of  the  "invariance"  of  these  functions  is  different 
from  that  of  the  second  order  invariants,  in  that  some  of  the  functions 
are  true  invariants,  while  others  can  only  be  considered  as  quasi -invariants. 
Four  methods  of  cascading  the  equations  are  discussed,  each  of  which  requires 
severe  restrictions  ©n  the  coefficients  a,  b, ...,g.  The  class  of  equations, 


for  which  each  method  will  result  in  a  termination  of  the  chain  after  a 
finite  number  of  ite rat ions ,  is  explicitly  pointed  out. 

The  final  extension  is  a  generalization  of  this  third  order  exten- 
sion to  the  n   order  linear  hyperbolic  equation  with  n  independent 
variables,  of  the  form 


(8) 


n 
U*l*2"-*n     f  !?!**%•  ••Xl-Ixl+1"  'x 


n 


+  f     a     .,  +  •  •  •  +  biux.  +  cu=  O 

i%i *k i\-"xHxi+r"xMxj+r-  ,xn 

Introducing  the  linear  operator    D,     so  that     (8)     becomes 

D(u)     =  0, 
the  number  of  identities  which  may  lead  to  a  decomposition  of      (8)     into 
a  system  of  two  equations  of  lesser  order  is  computed.     The  exact  number 
of  true  invariants  is  determined  to  be     n(n-l),  while  upper  and  lower 
bounds  on  the   number  of  quasi-invariants  are  formulated.     Theorem  V  proves 
that  an  invariant     h     is  a  true     invariant,   i.e.  invariant  under  the  change 
of  coordinates     u  =     /\(x^,X2*<> .  -x^  u  >  if 

h  =    3  %     +     atf0)    -     a..  .} 

for  some     i,j  =  l,2,...,n,     i  4  jj  while  any  other  invariant  is,   in  general 
not  a  true  invariant,  hence  a  quasi-invariant© 

A  discussion  of  earlier  extensions  is  included. 
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SECTION  I 


INTRODUCTION 


The  cascade  method  was  devised  "by  Pierre  Simon  Laplace  ,  and 

investigated  in  complete  detail  by  Gaston  Darboux.2  Earlier  ex- 
's In- 
tensions were  made  by  Uo  Dini<,J  and  J.  LeRoux  f   and  these  exten- 
sions are  discussed  "briefly  in  SECTION  VI  of  this  paper.  The  method 
itself  deals  with  the  linear  hyperbolic  equation  with  variable  coef = 
f icients  which  has  the  form 


M  X(u)  -    ^u   *  a(x,y)  ^u  +  b(x,y)  ^u_+  c(x,y)u  =  0 

<)xdy  5x         z7 


where  the  coefficients  ap  b,  and  c  are  analytic  in  a  certain  domain 
D,  or  at  least  as  many  times  continuously  different iable  as  we  feel  is 
necessary ,  and  u  —  u(x,y)  is  a  real  function  of  real  variables. 

Laplace's  cascade  method  begins  with  the  introduction  of  what 
are  now  called  the  two  Darboux  invariants,  h  and  k,  which  are  defined 
by  the  relations 

h  -  ^a  ±   ab  -  c  , 
2x  J 

(2) 

k  =•  **  +•  ab  -  c  , 

Then  (l)  may  be  written  in  either  of  the  following  forms; 


(3) 


3   (du         N  > 

^r~     Jy"    +  au^  +    b  ( JT    +  au)  "  hu        =  °  > 

A.  (  **_  +  bu)  ♦     a(   ^JL  ♦  bu)  -  ku  =  0. 
07       «  ex 


We  may  then  consider  the  obvious  substitutions 


u,  au 

1     =  rr —     +     au  , 

d  y 


(h) 

u~i  =  i2J  b* 


and  further  reduce     (3)     to 


(5) 


_     +     bu.^  -    hu  =  0  | 


+     au         -     ku     ■  0 
1 


d7 

Now  if  we  should  be  so  fortunate  that  either  function     h     or     k    is 
identically  zero,  we  will  have  succeeded  in  reducing  our  original 
second  order  equation  to  a  system  of  two  first  order  equations, 
either 


1     =      jj  u       +     au  , 

(6a) 

±       +     bu       -  0  5 
^  x 


. 


or 


(6b) 


"X  *  }x 


bu 


^7 


au  ,  =0 
=1       J 


either  of  which  syatem  may  he  solved  by  quadratures c 

If,  however,  as  is  more  likely  the  case  neither  h  nor  k 
is  identically  zero, all  is  not  lost.  To  be  specific,  let  us  consider 
the  system  (6a),  as  the  details  are  quite  similar  whichever  system 
we  choose „ 5  Instead  of  (6a),  we  have 


u,  -  ^u 


au 


(6a') 


1 


f  k,   -  hu  =   0  . 


We  may  integrate  the  first  equation  of  (6a5)  to  obtain 


(7) 


u     —    e 


^  dy     +     X(x) 


where  X(x)   is  an  arbitrary  function  of  x<,     Substitution  of  this  ex- 
pression for     u     into  the  second  equation  of  (6a°)  yields 

^ady 


|il   +.    bu 

>  x     x        1 


he 


r        (ady 


X(x) 


=  0 


or 


ady 


(8) 


^, 


+     bu 


3>x 


ady 


«]_       dy     +     X(x) 


_l 


Taking  the  partial  derivative  of  both  sides  with  respect  to  the  vari- 
able    y     give  s 
ady 


+  i 


5» 

L.     +     bu 

^x  H 

J 


<$  log  h 
d  7 


*    b       *\      +  3b      u 

d*^y  ay  -J7      3 


■> 


+     bu 


jady 


u 


After  some   straightforward  algebraic  manipulations 9  this  becomes 
2 


a 

+      1 


1 

a  x 


(9) 


3  u 

] 


+     c. 


«!     =     0  , 


where     a       =     a     - 
1 


a 


log     h 


3y 


c_      =      c     - 


b,     and 


+  ^  b       -     b 

3"*  "Jy 


<?  log     h 

5  y        * 


We  observe  that     (9)     is  of  the  very  same  form  as     (1),  but  with  new 
coeff ients       a  ,     b-,  >  and     c,    .     Hence  we  may  repeat  our  original  pro- 
cess,  in  hopes  that  either  of  the  new  invariants     h,      or     k-|_     will 
perhaps  be  zero.     If  nots  we  may  proceed  to  evolve  a  new  equation  for 


the  variable       u  s     and  in  fact  we  may  iterate  the  entire  procedure  as 


5 


often  as  necessary  to  produce  a  chain,  or  cascade  of  equations,  in  the 
hope  that  at  some  point  the  iteration  will  stop  because  one  of  the  in- 
variants will  become  zeroc  This  will  enable  us  to  solve  a  first  order 
system  by  quadratures,  and  by  tracing  our  way  back  through  the  chain, 
we  can  easily  solve  the  original  equation  (l). 

As  indicated  previously,  we  may  do  an  analogous  procedure  with 
system  (6b),  producing  a  cascade  of  equations  in  the  "opposite  direc- 
tion". Darboux  poimted  out  that  if  we  followed  the  h  substitution 
with  a  k  substitution,  we  would  not  produce  a  new  chain,  but  would 
indeed  revert  to  the  original  equation  (l).  In  fact,  if  we  denote  (l) 
by  E,  denote  the  equations  obtained  from  E  by  use  of  the  functions 
h,  h^,  hg*  ooo  by  E,,  Eg>  E~  00o,  and  denote  those  obtained  from  E  by 
use  of  the  functions  k,  k±f   kg,  „00  by  E  _x,  E_g,  E,,  ...  .,  our  chain 
of  equations  appears  as 

ooo,  E  „E    E„E.E 

If  we  should  take  any  E   obtained  by  use  of  an  h  ,  function,  and 

n  2i-l 

attempt  to  obtain  a  new  equation  using  a  k   function,  we  would  in  fact 

produce  equation  E  ,  0  The  same  relationships  hold  for  the  E   equations . 

n-±  .n 

Darboux  discusses  these  and  many  others  points  regarding  the 
Laplace  cascade  method,  including  the  nature  of  the  invariance  of  the  h 
and  k  functions,  periodicity  of  the  h  and  k  functions,  and  the  form 
of  the  most  general  solution  obtainable  if  the  chain  terminates  after  a 
finite  number  of  iterations  in  either  direetion„  It  is  the  purpose  of  -X 
this  thesis  to  extend  this  cascade  method  to  larger  classes  of  equations 
and  to  carry  out  similar  investigations  regarding  these  new  applications 
of  the  method. 


In  particular  we  will  first  show  how  the  cascade  method  can  be 
applied  to  systems  of  n  linear  second  order  hyperbolic  equations  in 
two  independent  variables,  with  n  dependent  variables,  discussing  the 
invariant  nature  of  our  substitution  functions,  the  general  form  of  a 
solution  when  the  chain  terminates  after  a  finite  number  of  iterations, 
and  the  solution  of  some  typical  systems «  We  will  consider  the  system 


(10) 


n 

r 


10 


au 


n 


a^ 


j-i 


13 


5  7 


n 


3   =  1 


c .  .  u. 
iJ   J 


0,  i  =  1,2,  e.  ,,n  , 


where     a^     =  *±.   (x>y)J   b_  =     b^foy),   c^     =   c±j(x,y)     are  continu- 
ously differentiable  as  often  as  necessary,  and  the     u^  =  Uj_   (x,y) 
are   real  functions  of  real  variables,,       We  will  first   express  this  equation 
in  matrix  form,  letting     A  =     (aij)*     B  =  Ov?)*  C  =    (c^j)   » 


U  = 


Un 


U, 


(o)  = 


01: 

0 


0 


0  = 


"  00 

•  • 

.    0 

00 

•   c 

•   0 

«  • 

• 

00     , 

•  • 

0 

__ 

Then  (10)  may  be  written  as 


(ID 


U 


xy 


AU 


BU    +  CU  =  (0) 


where  U   = 


&      3ft 
0  X^ 


~l 


ax 


U 


=  [' 


) 


u 
n 


5 


a  . 


Introducing  the  two  Barboux  n  x  n  matrix  invariants 

H  ■■  A   -h  BA.  -  C  , 


K  — 


4- 


we  will  be  able  to  proceed  with  development  of  chains  and  solutions 
by  the  cascade  method,, 

Our  second  extension  of  this 
order  linear  hyperbolic  equation  in  thr 
the  form 


will  he  to  a  single  third- 

indepeadent  variables,  of 


(12)    uxys 


+■  a(x,y,  z   )u   + 


b(xj,y,z  )u    +  e(x,y,  z   )u 

•*-Z  •*»/ 


-I-  d(x,y,z  )u   +  ©(x,y,z  )u   f  t(x,y,z   )u  +-  g(x,y,  z  )  u  =. 
*  »  z 


Since  we  will  he  dealing  in  three 
that  three  chains  will  result,  in 
and  the  wz -direction**.  Also,  sine© 
will  see  that  eighteen  "invariant" 
will  note  that  the  invariane©  ©: 
nature  as  that  of  the  corresponding 


=.  0 
independent  variables,  we  will  see 
Mx -direct  ion",  the  "y-direction", 

the  equation  is  of  third  order,  we 
functions  must  be  introduced „  We 

functions  will  not  be  of  the  same 
functions  in  the  second  order,  two 


must  be  placed  on  the   coefficients     a(x*y>z),    „    <>    »  9  g(x5y3x)  before 
a  chain  of  equations  can  be  developed* 

Finally  we  will  generalize   the  results  indicated  in  the  preced- 
ing paragraph  to  the  single   linear  hyperbolic  equation  of     n         order  in 
n     variables*  which  has   the  foim 

\ n  n  i 

_L* _     ♦     2:      ai    *        s 


n  ..    n~2 


+ 


(13) 


a„  . 


a 


1J  ^  u 


i5j  =  1        a^...^^  axL+1...axj_1dxJ+1#..dxn 

n  _ 

b,  ^ 

+     cu  =  0   „ 


«       o      o 


+  Z    bi    a  u 

In  this  section  we  will  introduce  an  operator  notation  which  will 

th 
greatly  simplify  the  calculations  for  n    order  equations  *  This 

will  enable  us  to  predict  the  number  and  form  of  the  identities s   and 
the  number  and  form  of  the  invariants  corresponding  to  each  identity 
related  to  the  n    order  equation o  From  this  we  will  be  able  to  in- 
dicate the  conditions  necessary  to  produce  a  cascade  of  equations* 
which  will  enable  us  to  tell  when  a  reduction  in  order  will  be  possible „ 
In  our  concluding  section  we  will  discuss  briefly  the  extensions 
made  by  Darboux*  Dini*  and  LeRoux  to  systems  of  second  order  equations 
in  one  dependent  variable  <,  to  a  single  second  order  equation  in  n 
indepent  variables ,  and  to  a  single  nth  order  equation  in  two  ^dependent 


variables »     These  extensions  are  treated  more  extensively  in  the  works 
of  these  men  noted  in  the  bibliography. 


SECTION  II 

AH  EXTENSION  OF  THE  PRODUCT  TJSTEGBAL  TO 
NON -HOMOGENEOUS  SYSTEMS  OF  THE  FIRST  ORDEE 

In  Section  III  w©  shall  show  how  the  cascade  method  can  he  applied 
to  a  system  of  second  order  equations.  The  end  result  we  hope  to  obtain 
is  a  reduction  to  two  first  order  systems, one  of  which  win  be  homogeneous, 

while  the  other  is  non -homogeneous .  The  solution  of  the  homogeneous 

6 

system  can  be  found  by  employing  Volterra's  product  integral  .  As  yet, 

however,  this  concept  appears  not  to  have  been  extended  to  non -homogeneous 
systems •  *  It  will  be  necessary  for  us  to  do  this  now,  to  enable  us  to 
solve  completely  the  second  order  systems. 
Consider  the  system  of  equations 

n 

^  =  y 

ax       j~-   aJi  uj  T  f±j 


(1) 


u^b)  =  uio  i  ss.  1,2,  »..,n, 


where  the  a   and  the  f,  are  given  continuous,  single -valued,  bounded 
functions  of  the  real  variable  x,  on  some  non-empty  interval  of  the  reel 
line,  b  ^.  x  =i  c,  and  the  u.   are  n  given  constants.  First  let  us 
write  (l)  in  the  notation  of  matrices.  We  denote  the  row  matrices 


UCxJ^Cu-l  Bg,,, 
F(x)— (f^  f p  •  .  •  fn) 

Uo  *  (U10  ^20  .  .  .Uno> 


10 


11 


and  the  square  matrix 


A  =    (aij) 


n  x  n  ' 


It  will  be  necessary  to  assume  that  the  matrix  A  is  non-singular, 
that  is  that  the  equations  of  system  (l)  are  linearly  independent. 
In  the  notation  above,  system  (l)  can  he  written  as 


(2) 


dU(x)    — 
dx 

U(h)  —  U 


U(x)A(x)  +■  F(x)  3 


o  • 


Let  P   he  any  partition  of  the  interval  [/b,ej9   such  that 
*  ~  x0  <  x^  <  ,  0  „  <  3^  =s  *,  and  let  j¥   he  any  point  in  the 
interval  /x^.      x^f.     We  then  define  U.,   hy  the  relation 


U. 


(3) 


U  ^  A(  \)   (x,  -  x^  )  +  U^F(  ^)(xK  -  XK4) 
U^AC^Kx,  <_a)+  i]  +  F^Ux^-x^). 


Here  I  is  the  n  x  n  identity  matrix  with  ones  on  the  diagonal  and 
zeros  elsewhere o  From  the  form  of  this  recurrence  relation  (3)  we 
readily  obtain  that 


w 


•n 


j«^+l 


+   F(<fm  )(xm  -  xm.x  )  f 


12 

Hext  we  consider  a  sequence  of  such  partitions,  ]!*.[  such  that 
as  m — »  oo  ,  A  x  -  xK  -  x,,  .  " — *"  0„  Sine©  all  the  functions  con- 
cerned are  continuous,  we  may  proceed  to  the  limits 


(5)      U(c)  -■   lim   U^  =  Uc  £  A  4    ]f(^)Ja  d^, 

h       \>  j 

e      c 

In  this  expression,     A  =   UA(^  )dJ)f  I  }   is  the  "right" 

h      h 

c     Tc 
integral  of  Tolterra,  while     C  2T($)  T  A$^   is  the  row  matrix  of 

term. "by  term  Riemann  integration  of  the  elements  of  the  row  matrix, 


F(^)   [  A. 


In  his  discussion  of  product  integration,  Schlesinger  proTed  the 

8 

following  identity s 


I 


q  *- 


T  <1 

C     ,      for       p  <.    s   <   q. 

s 


Using  this  we  see  that 

.1 


l   *   =  J  A   J   A   , 


h  <  ^  <  c, 


h    ^ 
and  hence 


!,-(C«)"I>. 


(6) 

Putting  (6)  into  (5)  we  obtain 

(7) 


1c        c       I  J  \-l  c 
A  +   C  F(^)/  |a\  Ta  dJ  ^ 

h       h      \b /  *> 
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r 


Since     A  is  constant  with  respect  to  the  variable  of  integration,  and 
"b 

since  the  scalar   &  "f  commutes  with  every  square  matrix,  we  may  factor 
this  product  integral  to  the  right,  and  write  (7)  as 

/  3  r1  Lc 


(8) 


U(c)  = 


u° +  5  r(1?)(I  AW 


■""Id 


This  function  U  which  we  have  derived  is  a  function  of  the  end  point, 
c,  of  the  interval  /"h^c/o  If  we  vary  this  end  point  in  any  interval  in 
which  the  functions  a. .  and  f*  remain  continuous,  single -valued  and  ■ 
hounded,  then  U(c)  becomes  a  function  of  a  real  variable  which  we  may 
call  x,  and  hence 

u(x)  -   juo+  ^  *{))  {  \    A\d^ 


(9) 


I 


x 


We  assert  that  (9)  is  the  solution  of  (2)  and  hence,  element-wise,  the 
solution  of  the  non-homogenous  system  (l).  To  see  this,  we  observe  first 

T  "^       tm        9  *  /   X  \-l 

that  when  x  —  b,     A  -  "£   ,  while  \    F(^)/l  A  £  "C  =   0. 

b  h      \h    ) 

Thus  U(b)  "=-  U  .  Next  we  note  that  each  row  of  the  matrix   I  A  is  a 

,,.     .  .  .  .  10  b 

solution  of  the  matrix  equation    d  »  s.  U(x)A(x)   .  Therefore 


Thus  differentiation  of  (9)  yields 


Ik 


J  A\   J  A   +   tJ0  Jf(| 


•1  Vx 

A4« 


A  \A  = 


=  F(x)  +  U(x)A(x), 

which  proves  that  (9)  is  the  solution  of  (2). 

11 
Schlesinger  showed    that  if  A  is  a  matrix  of  constant 

elements,  then 


a   —    e  ■ 


fc 


Therefore  in  the  special  case  when  (l)  is  a  system  with  constant 
coefficients,  (9)  takes  the  more  familiar  form 


(10) 


U(x)  = 


and  this '  solution  can  be  Terified  even  more  readily  by  direct  dif- 
ferentiation. We  note  the  analogy  "between  (12)  and  the  solution  of 
the  single  non-homogeneous  first  order  equation 
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g  ^    a(x)  u(x)  ^-  f(x) 


«(*)— B  uQ 


This  equation  has  the  familiar  solution 

rX 


u(x)  =     e  "         u.  t     ,    Jb     l 


0° <  I 


"b 
a  form  which  we  could  obtain  from  (12)  hy  talcing  the  transpose. 


SECTION  III 

SYSTEMS  OF  SECOND  ORDER  LINEAR 
HYPERBOLIC  EQUATIONS 

A.  Consider  the  following  system  of  second  order  linear  hyperbolic 
equations  j 


(1)    d  \ 


a,. 


Ju- 


n 


J*l 


u 


+  x 


au' 


n 


+  Z*Vj  *  °J 


J«l 


1,2,...  n^ 


where  the  coefficients  a*-,  "b^j,  and  Cj*  are  all  real  functions  of 
the  real  variable s  x  and  y,  continuously  different iable  in  both 

variables  as  often  as  necessary,  and  u^s  u^(x,y)  are  real  functions 

22 
of  the  real  variables  x  and  y.  To  put  this  system  in  matrix  notation  , 


let  A-(a±1),  B  -(biO,  C-Cc,.),  U« 


LU 


ij 


"la- 


ws obtain 


uii    (o)  •  r° 

u^  '  • 


% 


and 


(2)      U    +   AU   +"   BU   +■   CU 
xy        x        y 


(0) 
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where  U__„  =  52U  _, 


r^i1 
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dx^y 


^  x  5  y 


,  etc.  We  will  assume  that  the 


matrices  A,  B,  and  C  are  all  non-singular. 

We  introduce  the  two  Darboux  matrix  invariants 


(3) 


H  =*  A    +■    BA   «   C  . 
x  J 


K  m  B    +.  •  AB    -  C  . 


Then  we  note  that  (2)  can  he  written  in  either  of  the  forms 


HU    -    (0)j 


CO 


^pU        j.      Au\    +•      Bf   Ju     -h     Au"\      - 

TZVTT  J  \J7  } 

__,  (  %3L      +-       ^    +      A  /  ^  U     -h     Bu)      -    KU     —    (0) , 


If  we  consider  the  substitutions 


Ul  =   U7   +   OT  , 


U-l  = 


"*   + 


BU 


then  {k)   "becomes 


(5) 


»1  + 

X 

BU1 

-       HU      -     (0); 

u-l  "*■ 

AU-1 

»       KU      «■    (0)  ( 

H     =  0 

or    K 

—   0,  where     0  = 

O  0 

o  o 


our 


o  o   o 


n  x  n 
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original  system  (2)  will  be  reduced  to  two  systems  of  first  order 
equations,  which  can  be  solved  by  the  product  integrals  of  Section  II. 
We  will  make  the  assumption  that  none  of  our  coefficient  matrices  are 
singular,  except  the  matrix  0. 

Suppose  that  H  EL  0  .  Then  we  have  the  following  systems  to 
solve  ; 


(6) 


IL 


U 


BU. 


U, 


=   (o)  , 


+■ 


AU 


i      y 

The  first  equation  of  (6)  has  the  solution 

x        .-1 


(7) 


ui   = 


B(  If  ,y)df   T  (y) 


where   Y°  (y)  is  a  column  matrix  of  arbitrary  functions  of  y,  and 

B(  ^  >y)     d-^F  )    is  "the  inverse  of  the  Volterra  product  integral 
of  B.  Knowing  U  }     we  may  then  proceed  to  integrate  the  second  equation 
of  (6),   to  obtain  the  solution 


(8) 


U(x,y)  =  (  A(x,  71    )d7?j    X°  (x)   + 

'  ?       It  \~l         1 


where  X°  (x)  is  a  column  matrix  of  arbitrary  functions  of  x. 
If  H  3=   0,  but  K  =  0,  we  can  solve  the  system 


(9) 


U 


-1 


y 


AU.X  = 


U 


u. 


+    BU 
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by  the  same  methods  utilized  to  solve  system  (6),  to  obtain  the 
solution 


(10)  U(X,y)«     (f  B(f  ,y)d^J  * 


*°(y)  + 


(#  lf+I   []>(!^)M    (f  A«rt)^  *°(5>   4  J 


Illustrative  example    No, 

Consider  the  system  of  three  equations  in  three  unknowns 

"5^1  +.  7  o^i  +  y«y  5^  +y  ^*3  x  ^^i-J yai +  y®7^^    +     ^3s0j 
2x5T        5~*         JaT"      ^x      Jy~ 

(a'J  U2-e7^u1  +  y   ^  ^  ^  Ju^  ^     -  6^+ yUg  +  e^^  —  0 
^xj"y  J"x  £T"         Jx~   Jy~ 

^^^^l+^^+y®7^ +  ^u3+y2ui+°7u2+ye7u3 -°    • 

dxjy  ^*"~    ^~    V*   ^y 
This  system  can  he  written  in  the  form  (2)  with  matrix  coefficients 


A  = 


-eJ 


ye" 


y 


Q 


y 


I 


^ 


ye 


1 


B  ^sTl  o  o 

o  1  o 


o  o 


A 


If  we  compute  the  invariant  H,  we  find  H  —  0.  Thus  we  may  reduce 
to  the  two  first  order  systems 


w 


¥x   +   \    «  (0)  } 

x 


<V  + 


AU  ■  U 


1  • 


The  first  system  of  (b)  has  the  immediate  integral 


%  = 


-x 


e      Y(y) 

where  Y  (y)  is  a  column  matrix  of  three  arbitrary  functions  of  j. 
Using  this  expression  for  U-,  in  the  second  system  of  (b),  we  obtain 
the  solution 
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l(*ri  s  f  Ta(?  ty)  *  J  X(x)4UA(^)d^  e  Y(^)d^ 


where  X(x)  Is  a  column  matrix  of  arbitrary  functions  of  x< 
B»  Consider  now  successively  the  substitutions 

Us    f\   w'j 

«■  f  (y») 
x  -  y' 


C 


where  ^^  is  an  n  x  n  matrix  of  twice  continuously  differentiahle 
functions  of  x  and  y,  such  that  |  A|t^  ^    9   and  U  « 


u£  =  ui  (x^y).  Under  the  change  of  variables  U  ®  Au  ,  (2)  becomes, 
/yu^-h  (AA  +  ^y)U^  +  (B/\  *  Ax)Uy  + 

+  (CA  f  AAt  4-  B  Ay  +  ^xy)  u'  =  (0). 


Multiplying  from  the  left  by      A°      we  obtain 


U 


xy 


-1 


-li 


+   (AaA^^  Ay)Tj<    4-  ( AT bA  +•  A  /\  )u«     4- 


x'    y 


+  (  A"°T!  A  +  A-1a  Ax +■  A"^Ay  +  A"1 /\  ^u' ^  (o), 

.  is  of  the  same  form  as     (g)s 


(11)      U'        -4-      A'U'X     +    B'Uy     +       C'U'    ~   (0); 
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where  A*    =    A*1  A/\+  A"1  Ay 

(12)  b'  -   A^BA-t-A^Ax^ 


c*  -  A1  c  A  +  A"1  A  ax  *  A3B  Ay  *  A'A^ 


Let  us  compute  the  values  of  the  corresponding  Darboux  invariants 
h'  and  K1  for  (11). 


_'    —         •  •    •        » 

H     ~      Ax  +  B  A  "  C    m 


=  (a^aa*  A"\r)  +  (A^ba+a^a  kA^aa  +  a^j)- 


X  x 

-U  a     .      A  -1.    a      .    A-l  _  a  .      A  -1 


-(  A    c  A  +  A"^  A  +  A1  b  A     +  A   A  ;,  ) . 

x  y  xy 

Consider  .  -l 

AA    =  ij 
(13)      >*  >*         >*  ^~      " 

^A1^  ■    a"1    MA"' 

Using  this  identity1^  we  may  differentiate  the  first  term  of  H 
to  ohtain 

-1  k  -1      a         a  -1  i  "fa      a-1 


h  *  -AN    A   aA+    A  a  A*  A"  a  a    -AWa,  ♦ 

X  X  X  i 

t  A'V^y  -f-  A_1BA  A  +  A1  B  Ay  *  A"1  A,  A"2A  A  +  A*  A* A  Ay 
-AaAxy  -  A"XC  A  -  AlA  Ax  -  A"1  B  Av  ~  A'1  A 


?Asy  -  AXC  A  -  A'lA  Ax  -  A"1  B  Ay  -  A"1/  \ xy 

=  A"x(ax+ba-c)A 

=  A"1  HA. 


In  similar  manner  we  may  compute  the  value  of 
I*  «.  B*   +  A'  b'  -c°  =  A^kA. 
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We  see  then  that  this  change  of  variables  U  »  /\  Tj'  merely  produces 
a  similarity  transformation  on  the  matrices  H  and  K,  and  indeed,  if 
/\     is  chosen  so  that    /\  H-  K  A  t   then  h'  «  H,  and  should 
A  K  =  K/\    ,   then  K*  ~  K. 

If  we  make  the  change  of  coordinates  x  ^TW,  y  r  (f  (y*), 
then  the  resulting  equation 

Ux'y'  +"  V   (^)AUX«  *  <?  (xOBUy.  +  <p'(j!)  ^(yOCU  =  (0) 
also  has  the  form  of  (  2  )  •  Here 

a'  -  qAyO  a 
b'  -  cf)  (x)  B 

and  the  corresponding  matrix  invariants  are 

k'  -    (?'  (sj   ip'  (yj  K, 

Finally  the  change  of  coordinates  x  -  y  ,  y=  x   merely  has 

i       » 

the  effect  of  interchanging  the  invariants,  so  that  H  **  K,  K  s  H. 

Thus  we  see  that  the  term  "invariant"  is  correctly  chosen,, 

C.  In  the  likely  event  that  neither  H  nor  K  is  0,  let  us  consider  .> 

the  system 

Ul  =  uy  +  AJ3J 


(lh) 


U^^-  BUX  =  HU 


Solving  the  second  system  of  (1*0  for  U,  we  oh  tain 
(15)  ,  U  «  H"1  U^  •»-  H"1^  ( 

Utilizing  identity (13)  we  further  obtain 
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(16)  Uy~  -H"1  HyH"1^  4  H"1  U^     -    H"1  H^B*"1^  +  H'T^  Vj+S'S&w  * 

Then  mutliplying  the  first  system  of     (ik)     from  the  left  "by    H,     and 
substituting     (15)     and     (16),  collecting  terms  we  have 

U.         -t-       A^    +      B^        +     ClUj_    =     (0), 

where 

^      =       (HAH*      -    HyH"     ), 

(17)  »1     -    »  j 

Ci    *»    (EA*rl  -  HyH"1  )B    4-  By  -  8  =    AjBj.  +-B^  -  S  , 

This  equation  is  again  of  the  form  of    (  2 )  and  we  may  iterate  our 
process  in  an  attempt  to  reduce  this  system  to  two  systems  of  first  order*. 
Computing  the  matrix  invariants  for  system    (17)     we  have 

%    *    Alx  +  *X  Ai    -Ci 

»  (HAH"1     -  HyH"1)^.  +    B(HAH~       -  HyH*"     )  -  (HAH-1  -  HyH "fjB  - 

-  By  ♦  H 
S*8H-K«AX4AB-BA4-  B(HAH~1  ^*yaf     )  -  (HAH"  -  HyH      )B  + 

♦  (HAH"1),  -  (HyH"1)^ 

K-l*    H. 
In  a  similar  manner  we  may  solve  the  system 


2k 
to  derive  the  system 

(19)         U-1X7   *     A-1U-1^  +  B-1U-1    +  C-1U-1     ~     (  °  * 


where  now 


A.!    *     A 


B^    =    KBK"1    -    K^K"1 

C.x    —  (KBK-1  -  Kjt^)A  -  K  +    Aj.  = 

*B-lA.!+    A„^     _K^ 

For     (3  9)     our  matrix  invariants  will  be 

-1  ^5 

K_x    =  2K  -  H  -  Byi-BA.  =AB+A(KBK"1  ~  K^"1)    -  (KBK"1-  K  K-1)Af 

If  we  designate  the  original  system  of  equations  (  2  )  by  E, 
we  have  "been  able  to  transform  E  into  a  system  of  similar  form,  which 
we  shall  designate  E^,  "by  the  change  of  variables  U,  —  U  +■  AU,  and 
into  another  such  system,  to  "be  designated  hy  E  , ,  "by  the  change  of 
variables  U^  »  Ux  •»-  BUo  Provided  that  none  of  the  resulting  H  or  E 
invariants  is  non-singular,  w©  a&y  continue  this  procedure  in  both 
directions,  establishing  a  chain  of  systems 


o   o   • 


C     •     O 


E-3>E-1>  E>  Ei>V 
The  question  naturally  arises,  if  after  making  the  substitution 
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TJ^  as  U_  +•  AU  and  di  scorer  lag  that  H  ^  0  f  t$q   transform  system 
(  2)  into  system  (1?),  what  results  if  we  then  make  the  substitution 
Ul«,l  ~  Ut.  "*"  Biui  *  »  Presumably  we  shall  again  obtain  a  system  of 
the  form  of  ( 2 ),  hut  what  exactly  are  the  coefficients  of  this  new 
system  in  relation  to  the  original  system  (  2  )  I    Utilizing  the  expres- 
sions for  the  coefficients  of  (19)  and  recalling  that  KL=1,  we 


Al-1 

*      Ax    -    HAH"1       -  HyH*1 

Bl-1 

—    HHH"      «,    Ej"1 

cl«l 

*  B1-1A1-1  +    A^     "    H  m 

+  BA  -  EH"1!     -  AH"1  h  Jh"1   4 
J  x 

+  H^h-^-1^  HyH^H  H"1  -  H     H"1  -  H. 


consider  the  substitution  XL  ,  =  Au  9  and  for  <\  let  us  take  H, 
I.e.  U-,  -  «  HU  „  Irom  (12)  we  see  the  resulting  system,  again  of  the 
form  (2),  has  the  coefficients 

A'  a    H-l  a         H     4-    E~hL    —     A 

J.-1  y  > 

B*  ■*    I-3**-     H    +  H"1!     =     B  , 
1-1  x  ) 

S'  —    H"1  61-aI     +    I"1  A1-1HX    +    H"1  B^^Iy  +  H'1    H^y  »  C# 

Thus  we  hare  transformed  into  a  system  which  is  exactly  equivalent  to  sys&em 
(  2),   since 
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I»   l'-    1*1^1, 


"1-1 
^-  hence  H  ~  H,   .. 


This  result  may  "be  immediately  generalized  as  follows;     If,   in  system 
EA    we  make  the  substitution    ^j±i  ^^ty^A^i*  we  Will  obtain  system 
E±+l,  but  if  we  make  the  substitution    Ui„i~  ^ix  +  Bi^i>  we  will  ohtain 

the  system    E     _,   for  all  i  —  ,,..,    -2,    -1,   0,   1,   2,...    .     The  foregoing 

i— x 

tacitly  assumes,  of  course,  that  the  corresponding  Hj_  and  K^  are  neither 
equal  to  0. 

Remark  #1;  We  may  define  an  exponential  matrix  ©  ,  where  X  is 
any  n  x  n  matrix,  "by  the  formula 

e  *  i  4-  z  +  j_  +■  jr_  +  ...  -  i  +  2-  _x_  , 

2     3J  i*l  ii 

Now  unless  the  matrix  X  commutes  with  its  x  -  derivative,  we  cannot 

X    X 

assert  that  d6  ■  Q.     dx  ,  since 

dx        dx 

d_C   =  I  4-i  (dX  X  t  XdX)  +.  1  (dXx2  „ax  X4- X2dx:)   ooo 
dx  2  ax     dx    3<  dx     dx      dx" 

But  if  X=-Ay,  say,  where  A  is  an  n*n  matrix  of  constants, 

and  y  Is  a  scalar  matrix,  then  dX  —  A  dy_  ,  hence 

dx     dx 

XdX  »  AyAdy^  =  Ad£  *  Ay  r  dX  x  and  thus  d_(  e  X)  =  eX52, 
dx     dx    dx      dx  '  dx        dx 

Then  consider  the  special  case  of  system  (6),  when  the  matrices  A  and 

B  are  "both  matrices  of  constants.  An  integrating  factor  for  the  first 

fedX    ,J3X    mu  44- 

system  of  (6)  would  be  c  u   =  C  •  Then  we  may  write 
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'  _    "1 
H  —  H  — -  H   H-,  ™H , 
1-1  ) 

Hl-1 

**—■■    hence     h^H, 


ll-l 


• 


This  result  may  "be  immediately  generalized  as  follows:  If,  in  system 

E^  we  make  the  substitution  U.+1— U^i-AjiUi,  we  will  obtain  system 

E.^,  but  if  we  make  the  substitution  Ui.i*^i:ic  +  Bi^i^  we  will  obtain 

the  system  E   ,  for  all  i  ■=  ...,  -2,  -1,  0,  1,  2,...  .  The  foregoing 

tacitly  assumes,  of  course,  that  the  corresponding  Hj_  and  K^  are  neither 

equal  to  0. 

Remark  #1:  We  may  define  an  exponential  matrix  &  ,  where  X  is 

any  n  %  n  matrix,  by  the  formula 

X  2     3  t~    i 

©  «  I  4-  I  4-  X  4-  X~  +  . . .  -  I  + 


2     31  1*1  II 

Wow  unless  the  matrix  X  commutes  with  its  x  -  derivative,  we  cannot 

X    X 

assert  that  d  C  *  6.  dX  ,  since 

ax       dx 

d_C      =  I  4-  £  (<|X  X  +■  XdX)  4.  1  (dX  z24  ^  ax  I  4-  X2c]x) 
ft*  2     dx  dx         3i    dx  dx  ax 


O  0  o 


But  if  X=-Ay,  say,  where  A  is  an  n*n  matrix  of  constants, 

and  y  is  a  scalar  matrix,  then  dX  —  A  dy_  ,  hence 

dx     dx 

XdX  *  AyAd£  =  Ad£  *  Ay  =  dX  x  and  thus  dj  e  X)  =  eXdX. 
dx     dx    dx      dx  dx         dx 

Then  consider  the  special  case  of  system  (6),  when  the  matrices  A  and 

B  are  both  matrices  of  constants.  An  integrating  factor  for  the  first 

\Bdx   By 
system  of  (6)  would  be  e^   -  C  *•  Then  we  may  write 
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Bx        Bx 

e    u,  i  e    bu,  =■  (o) 


"X 


By 


fe(e    V=  «», 


♦  » 


Bx 


where     Y(y)  is  a  column  matrix  of  arbitrary  functions  of  y« 

-Bx 


Substituting  TL  —  a 

l 


Y  (y)  into  the  second  equation  of  (6)  gives 


-Bx 
e        Y(y)    =  Uy   +■   AU  ( 

\Ady  Ay 

An  integrating  factor  for  this  system  is  then     S^        -  €. 

Integrating  this  system  gives  tS&e  solution 

Ay  -  Bx 


u-  e"^ 


-  e 


(   Ay  -  Bx  s 

jc  x(y)ay    +-  x(x)J 

(Ay  f  Bx)    f  fe  ^  Y(y)dy   f    X(x)]  ^ 


where  X(x)  is  a  column  matrix  of  arbitrary  functions  of  x* 

Now  consider  system  (Ik)    where  no  longer  is  E  s.  O  ;  and 
assume  once  more  that  A  and  B  are  constant  matrices.  We  may  inte- 
grate the  first  system  of  (ik)    in  the  above  manner  to  obtain 

-Ay 


u  —  c. 


I  je^  u^ay  +■  x(x)J  # 


^dy  +-  X(x) 

Substituting  this  expression  for  U  in  the  second  system  of  ( 1^  ) 
leads  to  the  equation 

U  +  BU1  =•  H  e°^   je^  T^dy  4  X(x)  # 

Av   -1 
Multiplying  from  the  left  by  £  *   H     ,  this  becomes 


Ay  -i 

en     (u1   4-  bu 


f  Ay 
i>  ^)e    ui 


dy  4- 


28 


If  we  now     differentiate  "both  sides  with  respect  to    y    we  obtain 
( e      ah"      -  e  y  H~~%H~  )  (Uv   +    31^)  + 

+e^  h*(\  +  *%  +  Vl>  =  eAy  ui  ■ 

Multiplying  from  the  left  by  H  e."^   ,  and   collecting  terms  we 
finally  obtain 

\/       Al      ^*         ^       Bl     \        +         Cl      Ul        "        (Q) 

where  the  coefficients     A]_,  B]_    ,   and    C^are  exactly  as  obtained 
for  (17)o 

The  curious  thing  about  this  second  method  for  starting  a 

chain  of  equations  is  that  even  if     A     and    B     do  not  commute  with 

lAdy 
their  derivative B,  we  may  still  use  the  integrating  factor  c*  and 

d.^        kdy 
operating  formally  as  above,   pretending  that  the  equation    -^  -  £        «A 

uy 

were  actually  true,  we  will  arrive  once  more  at  system  (l?)  with  the 
correct  coefficients  A^  ,  B-^  ,  and  C^  «  Thus  w©  obtain  a  valid 
equation  from  an  invalid  operation,  but  an  operation  which  appears 
valid  on  the  surface,  and  is  completely  analogous  to  the  operations  per- 
formed on  a  single  equation  of  this  type. 

Do  If  the  chain  has  l)een  continued  to  the  i  4-  1st  system  of  equations, 
we  can  readily  establish  the  following  recurrence  relationships  for  the 
matrix  invariants  i 

(20)    Hi+laB2Hi  ~Ki  ~Aix+AiBi  -B^  +  B^AiHl  -  E^"1)  - 

-telAA"1  -=!_%%! +(H1A1H1"1)3C  -  (EU"1)^ 


Kj  ,  1   H. 
i  +  1    i  ) 


o  0  "*'-}      — -**J    9     "J        3  *  * 
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These  may  "be  solved  for  R.  and  K^     to  give 

(21)  *   1+lj  _x 

^^i+l       "     Hi+1       "     Al/     AiBi-BiAl+Si^l     ^i+l^i+l/^l  )        " 

"  (EIflAiKi+r     -  Ki+l/i+l"  )Bi+^i+lAiEi+l"  ^"(Ki+iyKifl"  )x 
From  these  we  immediately  ©Main  the  relations 

(22)  Hi^li'Ei-lS2Hi  "  Ah^  AiBi  "  BiAi*Bi(EiAiEi~  "  HiyHi"  }    " 


•(ViEi"     "  Hi/i"  )Bi  +^iAiHi"i}x  "  (HiyHi"\  j 


(23) 


Jy  j 

**»      ©  •  o      *-J    Jios   w%   JUo  C  y   m  a  o 

Finally  (  5 )  may  now  Tbe  written  more  generally  as 

^v+AU.1-E.(i=l)U.(l,1)> 

Considering  the  first  system  of  (  2k) ,   suppose  E  is  a  matrix 
of  constant  elements  c  Multiplying  "both  sides  "by  the  exponential  matrix 

Ti-1> 

B"1 

^-i^i-i      e       a  (e  "Ui). 

^x 

This  immediately  leads  to  an  expression  for  U  in  terms  of  U^  and  its 

derivatives  with  respect  t©  x^  namely 

TT    -1  ^=Bz  ^   ,^3x  -1  -%*    \     ,     Bx  -1  -Bx  n     n,  Bx  u  v 
d  s  ^x  ^~x    "Ji" 


B«  .  x 

<3  K  (see  remark  1,,)  and  solving  for  U,  ,,  we  obtain 

-1    -Bx  "\    3x 

J2_(e 
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fedx  fedx 

In  the  case  where  the  matrix    B     is  such  that     2L  (6U       )  =  B£ 

dx 

does  not  hold,  the  solution  for    U     in  terms  of    Ui    "becomes  much 
more  involved.     Solving  the  first  equation  of     (2^)     for    Li,  we 
have     n^  =  H^OJi     +BU±).     Then    U1-2  =  H^g"1^-;^    +    BQ^Js 

Hi- 2  (~L  ^Hi-l(uix+  BUi)  +■  ^i-l  *  (ui  *  BCT±d  •  Iterating  this 
process  we  see  that  eventually  we  will  obtain  U^jsU  in  terms  of 
U^     and  the     i     invariants     H,  B^, . • .H.   , .     Using  the  relation 


U±  =  ^ 


-ly       ^      Ai-l^-lj 


we  may,  in  a  similar  manner,  obtain  an  expression  for  Ui,  in  terms  of  U. 
E.  We  shall  now  give  consideration  to  what  may  occur  in  the  form  of  the 
invariants  as  we  iterate  along  the  chain  of  systems.  Perhaps  the  most 
natural  point  to  investigate  is  that  of  "periodic"  systems,  i.e.  systems 
such  that  after  J  iterations,  we  obtain  Ei  —  E.  We  shall  say  that  such 
systems  have  period  J.  Consider  first  a  system  of  period  1.  Then 
E-l=  E,  E1  =  h,  and  K-]_  =*»  K.  Hence  from  (20)  we  have  H  *  K,  and 

H  =-2H  -  E  -  A   4-AB-BA+  BCEAH"1  -  Eyfi"1 )  -(HAH"1   -  Hyl"1  )B+ 

+  im\  -  ■■(ByH"1)x. 

In  the  special  case  where  AB  —  BA  and  bCEAE"1  -  HyH^J-KHAH-i-HyH^jB, 
this  reduces  to 

(HAH"1)^  -  (HyH"1)x  -  ^  *   0 

which  has  the  immediate  integral 
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(26)         HAH"1  -HyH"1  -  A*-Y(y)^ 

where     y(y)  is  an  n  x  n  matrix  of  arbitrary  functions  of  y. 
If  in  addition  HA  ~  AH,  (26)  is  further  reduced  to' 

HyH'1  «    l(y) 
Hy    *  Y(y)  I  . 
This  may  he  further  integrated  hy  the  product  integral  to  give 

<27)  H  =  Q7  Yd^XU)^ 

where  X(x)  is  an  n  x  n  matrix  of  arbitrary  functions  of  x.  If  we 
select  for  Y(y)  a  matrix  of  constant  elements,  D  this  solution 
can  he  written 

(28)         h  =  ©*:  x(x)  ; 

Suppose  we  have  a  system  of  period  one  and  let  us  make  a  change 
of  variables  U  «■  f\  U  e  This  will  not  change  the  periodicity  or  the 
equality  of  the  invariants,  for,  as  we  have  seen,  if  H  *  K,  then 

H*  -  K1  hA*/^"  K  A  =  K1. 

Let  us  select  f\  however,  so  that  our  coefficient  A  ■■  0  ,  This 
simply  requires  that    Ay  _  =  -A  /\   ,  and  hence 

(29) 


A  -  J   {*®i    x  (*), 


* 


where  X  (x)  is  a  column  matrix  of  arbitrary  functions  of  x, 
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»  I  »•  '  t  I     '  t       *I         'l    t    '    I        I  *        I 

Then    E     =A       +BA       -    C  =    -  C    =    Z,     l?utS-B  >AB  -C  *S     -  C, 

a  y         y  * 

hence        B_  ■   0  .  If  we  make  further  assumptions  regarding  the 

character  of  A  and  B,  our  system  can  he  reduced  further.  Let  A  and 

B  he  matrices  of  constant  elements.  Then  (29)  can  he  expressed  as 

(30)  A  -  e     x  (x). 

Since  B  -  A  b  A  +  /\ "A*  »  we  ma^  a-^BO   attempt-  to  find  A  such 
that  B  »  0  •  This  gives 


(3D  f\   _  - 


e 


(y) 

Comparing  (30)  and  (31)  and  noting  that  e?eQ-  eP"'"Q  if  and  only  if 
PQ  "  QP,    we  may  take  for  our  arbitrary  matrices  X*  *  q     -   y  =.  ©"*v 
to  get  the  matrix 

(SB,  A-e"^*^ 

if  and  only  if  AB  *  BA. 

Using  (32)  to  change  coordinates  under  U  =  /\  U  ,  our  resulting  coef- 
ficients are  \. 


t 
A 

- 

0    J 

1 
B1 

•=? 

°j 

c' 

= 

< 

and  our  system  (ll)  has  the  reduced  form 


»       _«_» 


(33)  U    »  Hj 

xy 

which  Is  the  matrix  analog  to  the  telegraph  equations.  Since,  If  A  and 
B  are  constants  such  that  BA  »  AB,  it  follows  that  E  sss  Z,  we  have 
proved  the  following: 
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Theorem  I.  A  system  of  equations  of  the  form  (2)  haying  constant 


matrix  coefficients  A  and  B  can  "be  reduced  to 


form    U*    =■  iV 


"by  a  change  of  variables     U  =  AU       if  and  only  if  A3  =  BA., 
Illustrative  example  #2: 

-3x57  3x       ^x       dx       ^y~      3"7         ^  J 


}  u2  4-2^ul-  2*^2  +^u3     -  ^H-1^  +-12^u3  +  20ul  "  17u2  ^^u3  -  0^ 
(c)  dx^y   £x~     ^F"  7F"         Jy~"  <    ay         Jy~ 


S  u3+^>  +  ^u3  f  23ulf  2^u2  +  17^2.  f  6Ul  *  19U2^  79U3  "  °  ' 
^x^y^x       Jx       Jy      <3y  dy 

This  system  may  he  written  in  the  form  (2),  vith  the  constant 
matrix  coefficients 


A  ~ 


PL       3       5~ 

12-2       1 


0       1      k 


BsfY         2 


J 


28 


-2       11         12 


U 


17 


We  oh  serve  first  that  det  A  "=  -23,  det  B  «529,  and  det  C  =  -8,792,  so 
that  all  matrices  concerned  are  non-singular o  (Notice  that  B  »  A  and 
C  =.A^  -  I.)     Computing  the  values  of    H     and    K  we  ohtain 
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H  as  Av   4-  BA 


11    45  l49| 

20   -16   49 

6   19   80 


20 
6 


45  149 
-17   49 

19   79 


-   I 


K  =  Bv  +  AB  -  C   - 


45  149 1 
-16   49 
19   80 


10   45  149 

20  -17   49 

6   19   79 


0  0  0, 
0  o  ol  +■ 

0  0  0 


1  o  ol  *= 
0  10 

loo  1 


Thus  H  s  K  *  I,  and  we  further  observe  that  AB  *  BA.  The  conditions 
of  theorem  I  "being  satisfied,  (c)  reduces  to 

(d) 


l\ 


t  » 


xy 


h  u. 


Under  the  transformation  U  »  /\U  ,  however,  we  know  that 


-1 


k-l 


.-1 


h'  "=      A  lA=  ^  ^A-AA-t=H^  so  that  (d)  is  in  reality  a 
system  of  three  "uncoupled"  telet_vaph  equations, 


U*   ss  U* 
xy 


Each  equation  of  this  system  has  the  same  solution  u  ,  which  can 

"be  found  in  the  literature.  '   Thus  (e)  has  the  solution  U  =.  fu' 

u| 
u 

•   fa 

Having  determined  U  ,  we  must  then  compute  the  solution  U,  using 
the  relation  U  ■=.  /\  U  .  Equation  (32)  tells  us  that 


K         -(Ay  +  Bx} 


-y 


1 

3    5 

7     2     28 

2 

-2     1 

—    X 

-2    11    12 

0 

1    4 

2     2     17 
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hence  our  solution  is 

-(Ax  +   By) 

U   «    6  jjlm 

F.     We  consider  next  the  systems  of  period  two,    so  that     E2  —  E, 
and  hence     H2  «  H     and    Kg    =.    K.     Equations  (20)     and     (22)     tell 
us  that 

(3*0  2H  ^  2K  -  A^  -t-A^B!  -  B^    4-  tjjLXkjK*1  -  K,^"1)    - 

-(KA^"1  .  KyK"1)^  4-  (^K"1),  -  (KyK"1) 


X; 


and 

r"1    _  it  w-ll 


2K»2H-AX  +   AB-BA  +  B(HAH        -  HE"1)    - 


-  (BAH"1  -  HyB^jB  4-    (BAH"1^-  (HyH"1). 


Summing  these  two  equations  we  obtain 

(36)  r  -1-1  -1  -1       1  r  1  -1 

+jBi  miK      -  V    ^     "    ^l*      "  V    'Bi\  +     jB(HAH      -  Ijl  "V- 
-(HAE-1;  .  yf»>!  -  (yf)r-^)x  *  O. 

Looking  at  the  terms  within  each  pair  of  curly  "brackets,  we  note  that 
under  certain  obvious  conditions  of  commutativity  this  system  reduces  to 

(37)  ^'\^   (Ey^x  =  °  • 
Equation  (37)  bas  the  immediate  integral 

(38)  KyK"1  4-  HyE1   =   Y(y) 


36 


where    Y(y)  Is  a  square  matrix  of  functions  of  j. 

If  we  impose  another  condition  of  commutativity,  whereby 
EH  «  HK,   KyE  =■  BKy  then  an  integrating  factor  will  be  multiplica- 
tion  from  the  right  by  KH.  For  this  leads  to 


K^  +  H^K  -   Y(HK)  } 

J„  (HK)  =■  Y(HK)  . 
>7 


Hence 
(39) 


r^ 


HK- 


IJ- 


(Y(fc)dh>  2(x) 


where  X(x)  is  a  square  matrix  of  functions  of  x. 
These  functions  of  x  and  y  are  in  fact  determined  by  the  coeffi- 
cients A,  B,  and  C.  Thus  we  have  proved  Theorem  II: 

If  a  system  of  equations  of  the  form  (2)  has  period  two,  and 
the  following  conditions  of  currautativity  are  satisfied; 


**L  * 

A1K 

HA     as 

AH 

AB     as 

BA 

A1B1  = 

BlAl 

BlKyK"1 

as 

V 

■\ 

BHyH"1 

—        HyH 

-h 

KH  =   HK 


y 


then  the  product  of  the  two  matrix  invariants  H  and  K  has  the  form 
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HK  »   Y*  (y)  Z(x)) 

A  a 

where   Y  (y)  and  X  (x)  are  square  matrices  of  functions  of  y  and 
x  respectively  „ 

Since  H  is  non-singular,  there  exists  a  matrix  ©  defined  "by 
the  relation   H  =  @^  ,  so  that  <Q    =  log  H,    Uow  if  the  conditions 
of  theorem  II  are  satisfied;  and  in  addition,  (39)  takes  the  simplified 
form  HK  ~  I,  then  (3k)     can  be  "written 

2  H  =2S-  (lyr1^  , 
or 

r-1 


2  H  -  2  H"1  —  (H**H  ) 

y  x^ 


hence 


Furthermore,    if    ©    has  the  property  that      Q     ^Q     -       <?©  @)    , 
then    tt_     s   e       iO  =  H    <**&  ,   so  that 

Substituting     (kl)     into     (kO)     yields  the  equation 

£^y 

and  this  non-linear  system  of  equations  -will  then  have  the  solution 

(5)  —   log  Ho  Equation  (k2)     is  the  analog  of  Liouvilie's 

17 
equation. 
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G.     We  can©  now  to  a  consideration  of  the  general  form  of  the  solution 
when  some    Hj=0       H*  •£    0-    0  £j<i  .     Tliat  is,  we  it©rat©  th©  process 

•J       9 

forming  the  chain  of  equations  E,  E^  E2,  ...  E^  until  R±  s.  0  .  It 
will  then  he  impossible  to  form  equation  E.  , ,  hence  E.  is  the  first 
system  for  which  the  H  invariant  vanishes.  The  first  equation  of  (k) 
then  has  the  form 


Jfcta  +■  A±  uA-f  Bj  /  ^  +  A^A  =  (0)  . 


This  has  the  Immediate  first  integral 


fa  +  A^y  ff  x^y1 y^ 


3r 

where   Y  (y)  is  a  column,  matrix  of  functions  of  y  only.  Upon  inte- 
grating a  second  time,  we  obtain  the  solution 

vAx,j>{f  ^(x,  ^  )  d  )       |X(x)  f ;  '■;(:>) 

+J  fprf^)«),  (I  W*)     Y($>  J  45|. 

where     X(x)     is  a  column  matrix  of  functions  of     x  only. 
This  solution  is  of  the  form 

u^fxW     -f-     C7<B  YCyJdyj^ 

where  ^  and  (^  are  n  x  n  matrices  of  certain  well  defined  functions  of 
x  and  y. 

Using  the  iteration  scheme  of  paragraph  D,  we  see  that  the  solu- 
tion for  U  has  th©^  form 

+  a2(*«  +fve>  xdj)+  ...+^i  ( x(i)+  f^e  Ydy). 
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Here^jfl.,,  .oo5W.£  designate  certain  ¥©11  defined  n  x  n  matrices  of 
functions  of  x  and  y,  and  are  ecmposed  of  certain  products  of 
Hj   ,  Qa,   B,  and  their  derivatives.  We  note  that  the  column  matrix 
Y  always  appeal's  under  the  integral  sign,  in  the  most  general  ease. 
If  the  "boundary  conditions  are  such  that  Y  —  \P)  we  see  that  the 
solution  is  of  the  form 

and  this  is  the  most  general  solution  in  which  appear  arbitrary  functions 
of  x,  and  no  integral  signs. 

Conversely,  if  the  original  system  of  equations  has  a  particular 
solution  of  the  form  of  (45),  then  the  repeated  application  of  the  cas- 
cade method  will  certainly  lead,  after  a  number  of  iterations  not  greater 
than  I,  to  a  system  for  which  Eji  «•  O  .   We  will  now  prove  this,  i.e. 
if  (1*5)  is  a  particular  solution,  the  number  of  iterations  necessary 
■before  any  H-in^ariant  vanishes ,  is  not  greater  than  i.  If  we  substitute 

(45)  into  the  system  (.2),  we  will  obtain  an  expression  of  the  form 

From  the  arbitrary  character  of  X,  we  assert  that 

(46)  M   sO       J  =  0,1,2,. .„,  i*l  . 
If  we  compute  the  form  of  y^   ana  Jfc  .  -   we  ohtain 

d  -f      xl^x^y     ax     77    ^i; 


1 


a  y 


ko 


Equation     (k6)     implies  that 

w>  &l  +  A  a.  -  o 

c?  y  i  ' 

thus   *i  =    JO.^    +    A^.1+.     J     flfln.    AA^+Bf^^    A«A 

Hence  we  have  that 
(W)  ^i-1     *       A^wf        H*i# 

*y 

Recall  now  the  first  substitution 

<*9)  »i«  ^U  4-   Art 

If  we  substitute  the  expression  for  TJ  given  "by  (k5)     into  (^9), 

aid.  afcilize,  ($%)>  w®  se®  insaediately  that  the  resulting  expression  for 

U^  will  have  no  derivatives  of  order  greater  than  i  -  !„  In  addition, 

(i-l) 
(U)  informs  us  that  the  coefficient  of  X     in  the  expression  for  U|_ 

will  he  Hfl^  .  Whenever  H  vanishes,  therefore,  the  coefficient  of 
X^   '  in  this  expression  vanishes  also,,  That  is  to  say,  J.f  H  vanishes 
the  expression  for  U^  would  have  no  derivatives  of  X  of  order  greater 
than  i-2„  As  a  consequence,  repeated  application  may  lead  to  a  system 
E,,  for  J  <  i,  for  which  the  invariant  H.  s  O  j  otherwise  we  will 
obtain  an  equation  E^,  which  has  a  particular  solution  of  the  form 

u±  -  0.  x  . 

To  srammarize,  we  have  proved  the  theorem  III; 
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If  system  (S)  has  a  particular  Integral  of  the  form  (**•!?),  with 
Ut£  #  0   then  the  method  of  Laplace  will  lead,  after  a  number  of 
iterations  not  greater  than  i,  to  an  equation  which  is  integrable. 

Ho  Consider  now  expressions  of  the  form 

?  (1) 

0  ■  a  x *■  4  x  +  -.  1-  (xp    j 

•which  contain  a  matrix  of  arbitrary  functions  of  ac,  and  derivatives 
of  this  matrix  up  to  a  specified  order.  It  is  obviously  quite  often 
possible  to  express  U  in  terms  of  derivatives  of  an  order  greater  than 
specified.  That  is,  if  X  is  expressible  as  a  sum  of  matrices? 


where  (8^  (*»  •••  ^  ar®  certain  matrix  functions  of  x,  and  X^  is  a  new 
matrix  of  arbitrary  functions  of  x/  then  U  may  contain  derivatives 
of  Hi   up  to  order  1  +■  U   .  Conversely,  it  may  be  possible  to  reduce  the 
order  of  the  highest  derivative  appearing  in  U.  For  instance,  if 

u«  ^,(14  xl)  +  6(x-+ xl#)j 

then  the  substitution  Xx  —  X  +  X*  will  reduce  the  order  of  the 
highest  derivative  by  1. 

We  will  say  that  the  matrix  U  has  rank  14-1  with  respect  to 
x  if  the  highest  derivative  of  X  appearing  In  the  expression  for  U  is 
of  order  i,  and  it  is  impossible  to  reduce  this  order  by  substitutions 
of  the  type  described  above. 

We  assert,  that  if  the  system  (%)  is  the  first  for  which  the 
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invariant  H.  vanishes,  then  there  exists  a  particular  solution  of 
the  original  system,  composed  of  a  column  matrix  X  of  arbitrary  func- 
tions of  x,  and  derivatives  of  this  matrix  up  to  and  including  order  i. 
This  solution  is  irreducible  in  the  order  of  derivatives,  and  hence  has 
the  rank   i  +■  1. 

For  if  there  were  a  substitution  which  would  reduce  the  order 
of  the  highest  derivatives  to  i  -/J   ,  say,  then  there  would  be  a 
system  (Xj)  £  s  1  — jj   <  i,  for  which  this  new  expression  would  be  a 
particular  integral,,  But  this  would  imply  that  H,  —  0,  which  is  con- 
trary to  our  hypothesis.  As  a  result,  it  is  evident  that  if  (Ei)  is 
the  first  system  for  which  EL  5*  0,  then  no  other  system  (E.   ,  )  with 
positive  or  negative  indices  will  admit  a  solution  of  rank  i+1  with 
respect  to  x. 

The  results  of  the  preceding  discussion  apply  without  modifica- 
tion to  the  second  substitution,  i.e.  the  substitution  which  results  in  . 
the  chain  (E^)  (E  p)...  If  this  chain  eventually  terminates  at  some 
system  (E  ,)  for  which  K  .  =  0,  then  there  exists  a  particular  integral 
of  rank  j  +■  1  with  respect  to  y^ 


U*6Y  +  gjl  '+  ....+ej  1    (J), 


and  conversely. 


I .  We  are  now  in  a  position  to  construct  all  the  systems  of  a  given 
dimensionlwhere  by  the  dimension  is  meant  the  number  of  dependent  vari- 
ables) of  the  form  (2),  which  will  lead  to  a  general  integral  by  this 
extension  of  Laplace's  cascade  method.  Suppose,  for  example,  we  wish  to 
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construct  a  system  of  dimension  n,  which  will  Have  a  solution  of 
rank  i  4-1  with  respect  to  x»  Then  the  chain  must  terminate  after 
i  operations,  so  that  the  system  (E.^)  will  he  integrable.  First 
we  choose  arbitrarily  matrices  A.  and  B.  non-singular,  of  order  n. 
Then  the  equation. 

(50)  Hi   ■   A^  +  3±A±     -  d  »   (  O  ) 

will  determine  Cio  Then  a  particular  integral  is  given  by  (^3)» 
The  value  of  the  K^  invariant  for  (Ej^)  is  determined  by 

(5D  *1  *   Biy   -I-  AlBl  *  C±, 

and  then  the  relations  (17)  and  (21)  will  permit  us  to  calculate 

the  invariants  for  the  systems  (E^  _i)> .  «>  <>,  (E), 

Viz:      Hi.!  *  *u 

Bi-1  ~  Bij 

-1  -h 

Ai-1  -  Hi-1  (Ai  *  Hi-V  **  ^1^. 
(52)         .  h1d1  "I  k±   Hl-1  4.  h^'1  HjL^ 

Ki-1=  ai   "  Hi  "   Aix+  Ai^r-l  V?i^i.i 

+  "i^A"1   "  Ki    V1*  "  ^i-iV^V  Ki1)Bi-i  * 

y  y 

These  relations  also  permit  the  calculation  of  the  coefficients 
A,  B,  and  C.  From  the  expression  in  (52)  for  A^  and  B^,  we 
can  show  that 


hk 


-1 

A    *  (Hi.1Hi^..^51H)   'A^JI^..    HlH)  + 


i-1  i 

+  Z    (TT 

n=  1      \  jan 


B  =■   B 


b 


and,   having  determined    K, 

C   a   By  +■  A3  -  £  , 

Having  thus  determined  all  of  the  invariant  matrices  H  ,  0£j  <•   i, 

J 

we  may  then  compute  the  general  integral  for  U  "by  the  iteration 
process  of  paragraph  D. 

We  may  determine  in  a  similar  manner  all  systems  which  ter- 
minate in  "both  senses,  and  hence  admit  a  particular  solution  of  the 
form 

u  =   ^  x  +^xx'+  .•^1x(1ien  ...<^y(j) 

which  contains  no  integral  sign,  Th©  preceding  expression  has  rank 
i-*-l  with  respect  to  x  and  rank  J  -H  with  respect  to  y.  The  sum 
i  +  j  is  called  the  characteristic  number  of  the  equation.  This 
number  does  not  change  upon  successive  applications  of  the  Laplace 
Method.  In  passing  from  system  (B)  to  system  (E»  )  for  example, 
the  number  i  is  diminished  by  h,  but  the  number  j  is  increased 
by  the  same  amount;  the  sum  of  the  two  is  unchanged.  This  is  apparent 
if  we  consider  the  system  (E^),  with  invariant  1^*0.  This  is  of 
the  form 
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(53)     h{rr+ Al  u*)+ B* ( j-r  +  Ai  Bl)  =  (0); 

which  admits  a  solution  of  rank  1  with  respect  to  x,  and  of  rank 
i+j  +1  with  respect  to  y.  For  "brevity,  let  n  =  i  +"  J  and  con- 
sider the  substitution 


where  (Q)   *s  a  col^-1321  matrix  of  unknown  functions. 
Then     }tt  -  N  , 

o  *  3? 

hence    ^li-AjTFi  -  (Ja.47)"1  ^®  . 

S7  ^ 

Using  this     (53)     "becomes 

CO         £  { (I A^,"1    i®  J  ♦  B±  {  (J a^)"1    i|  >  (0)  . 

An  integrating  factor  for     (5*0   is  multiplication  from  the  left  by 
TB^dx.     Thus  we  may  write 

(55)  i       (TBidx     (lA^)"1      to)   =    (0)# 

Denoting    f^ta  •    (Udy)"    -    <*  ,      (55)  "becomes 


A   foe-1 4®}  =  (o> 
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which  may  "be  immediately  Integrated  to  give 

(56)  ^    .    Ot^     ,  ©sJcCY^    +     X1> 

where  Yi  >  designates  a  column  matrix  of  arbitrary  functions  of 
y,  and  XA  a  column  matrix  of  arbitrary  functions  of  x.  We  know, 
already, that   (56)     admits  a  solution  of  the  form, 

(57)  ©  =  Xx+  fi  Y    •»-     6^  +   ...finY(^ 

in  which  ©,  6i»»<>^  fin  are  certain  well-defined  matrices  of  functions 
of  x  and  y,  and  Y  is  a  colum  matrix  of  functions  of  y  only. 

If  we  substitute  (57)  in  (56)  and  take  an  arbitrary  numerical 
value  x,  we  find  that  Y    and   Yl  are  related  "by  an  expression  of 
the  form 


(58)         ,,  =AT  +  \x'  +  -+AnilT(ntl)J 

in  which  A  ,  A  ,  ...  A  -,  are   certain  square  matrices  of  func- 
tions of  y  only.  Substituting  the  expressions  for  <Q)     and   Y^  given 
"by  (57)  and  (58)  into  the  first  equation  of  (56)  implies  that 


r(n+l) 

In  order  that  (59)  he  a  valid  equation,  the  coefficients  of  like  order 


(59)  2fer+«iA...*ftl*Wf  =  <X(AY*AlYV..t  An41r(n+1> 

5y 
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derivatives  on  each  side  must  "be  equal.   [In  this  manner  "we  obtain 
the  system 


(60)  ^6=0C/\ 


4&  +  6  =  <*\ 

*  •  • 


(61) 


If  we  eliminate  the  |L  from  (60),  we  obtain  the  equation 


a  linear  equation  of  order  n  f  1  with  respect  to   OC  >  vhowr  matrix 
coefficients  contain  functions  of  y  only.  If  we  solve  (60)  for 
the$.,  we  obtain  the  values 

(62)  ®n~<XK<l, 

6  -  tfA,  -  1,  (oc  A  M  ...+  (-i)n  3  »  (ocA    A 

1       ^y  2  T-jn  \  n+y  j 

which  permits  the  determination  of  the  expression  for  ©  . 

The  relation  (58)  ,  between  the  matrices  Y  and  Y-  ,  per- 
mits an  arbitrary  choice  of  either  one  or  the  other  for  the  defini- 
tion of  the  value  of   ©  ,  by  either  (57)  or  the  second  equation 
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of  (56).  The  choice  of  the  former,  however,  will  give  the  most  general 
form  of  the  integral  for  Q)  . 

As  an  application  of  the  preceding  discussion,  which  is  com- 
plementary to  our  original  problem,  we  see  that  it  is  sufficient  to 
choose  n  +  2  square  matrices  of  arbitrary  functions  of  y,  namely 

}   1V"»    H-t-1   >   solve  the  resulting  linear  equation  of 
n  4-  1    order  in  one  independent  variable,  (6l),  for  the  matrix  OC  , 
in  order  to  determine  the  form  of  the  equation  (EjJ  .  The  formulae 

(62)  and  (57)  then  provide  the  general  integral  of  the  equation  (56) 
and  hence  (53)»  Repeated  application  of  the  substitutions  of  Laplace 
will  enable  us  to  determine  the  system  (K)  which  has  this  general 
integral. 

Now  consider  an  arbitrary  linear  equation  of  (n  +  lj    order  of 
the  form 

/  ^  (n±l) 

(63)  Mo(  4   M.QC  +  ...  +  M    tt       -   0  , 

1  n  4 1  D 

where   M,  M  . ..,  M    are  any  n+2  square  matrices  of  known  functions 

y  and  o£  .  is  a  square  matrix  ©f  functions  of  y  .  Let  A,  Aj^ ...  An+^ 

be  defined  by  the  identity 

a   '     Ait'       a      (n+-l) 
(ft)     *W  +  V  +  A^W  +-^An.lV 

-Mf  -  2(M1W)+  £?  (M2W)-...^(-l)n+1  ^(m  ,W>, 
where  W  is  any  arbitrary  square  matrix  of  functions  of  y.  (64) 
then  determines  the  A'3  in  ^orme  of  tne  M's>  for  ©sample 
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(65) 


iul 


In  light  of  equations  (60)  and  (65)  we  may  set 


n+1  > 


M 


11  + 


=   {-DJ 


+  1 


o{     e. 


n 


That  is,  if   OC  is  a  solution  of  (63),   then  we  may  determine  the 
functions   Aj   "by  the  identity-  (64),  and  the  Q      "by  the  system 
(60).  Upon  substitution  into  (57),  we  see  that  the  general  integral 
contains  the  column  matrix  Y   and  its  derivatives  up  to  order  n. 
Thus  we  have  proved  the  following ; 
Theorem  IY;  In  order  that  the  linear  equation 

M*'1  *&\  -  <»> 


for 


admit  a  solution  of  rank  n-4-1  with  respect  to  js   that 
general  solution  to  "be   of  the  form 


1         n   j 


such  that  the  solution  cannot  "be  written  in  analogous  form  in  which  there 
appear  fewer  derivatives  <gf  the  arbitrary  column  matrix  Y   ,  it  is 
necessary  and  sufficient  that  OC  s   considered  as  a  function  of  y ^satisfy 
a  linear  equation _qf_  order  n  4-  1  whose  coeff  ieJents  are  functions  of  ys 
and  that  OC  does  not  satisfy  a  similar  linear  equation  of  lesser  order. 


J.  The  results  given  in  the  section  have  "been  kept  in  close  analogy  to 
those  contained  in  Chapter  H,  volume  2  of  Leepns  Sur  La  Theorie  Generate 
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Pes  Surfaces  "by  Gaston  Darboux.  In  this  chapter  Darboux  accomplishes 
more  "with  the  single  second  order  linear  hyperbolic  equation  than  is 
available  to  us  with  a  system  of  such  equations,  due  to  such  causes 
as  the  non-eommutativity  of  the  ring  of  n  z  n  matrices  over  the  field 
of  functions  of  two  real  variables.  It  is  of  course  probable  that  some 
of  the  results  obtained  by  Darboux  which  were  not  considered  by  us  in 
this  section,  could  be  obtained  for  systems  "by  matrix  methods.  The 
reader  is  referred  to  this  very  fine  work  hj  Darboux  which  has  been  the 
inspiration  for  much  of  this  thesis . 


SBCTIOK  17 

HYPERBOLIC  EQUATIONS  OF  THUD  ORDER  IN 

teres  nm&Ewmm  variables 

Ao     Let  us  now  turn  our  attention  to  fc  near  equations  of  third 

order  in  three  independent  variables  1  have  the  form 

(1)  o^Cv)    =■    u^z     ^    an  -^    +    ct^.    +    du^  t-    eUy   4    fuz -t   gu  =    0. 

The  coefficients  a,  b,       .  f,  an  are  to  be  considered  as 

functions  of  2.,  aiiabl©  as  many  times  as 
wo  may  needc  the  manner  of  Laplace  dis- 
cussed previous  )  to  a  system  of  three 
first  order  equations.  .11  then  attempt  to  cascade 
the  equations j  in  the  hope  that  a;  be  number  of  iterations  the 
chain  will  terminate  .  invariants,  and  that  the  resulting 
system  can  then  be  reduced  as  aired* 

To  coMLence  this  operation,  we  must  consider  substitutions  of  three 
different  types,  n&Tcely 

(2)  u,   *      +    au  , 


(3)  -  ux  +   bu 

y       > 

-  u   +   bu   +■   cu   4-   du0 
-j.     yz      z      y 

Let  us  start  e  of  variables  (■"),  Using  the  definition  of 

the  operator  cJC  which  appears  in  (l),  w©  read:'     :ify  that 
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(5)    u    +-  "**!    4-  eux  +  dux  -  <A.  (u)  +*  k^  +  lx  uy  + 

jrz      z      y 


where  k.  s   a   +•   ae  -  e 

1       z  ) 

I   s  a   +   at  -  f,, 

(6)  ii       y  > 

*i  *  a   4-  «a  +  "ba  4*  a    -  g 

x     yz      J     z     yz 

are  defined  as  the  first  three  x -invariants.  It  is  apparent  from  the 
symmetry  of  the  equation  that  we  could  have  just  as  well  made  the  substi- 
tutions 

(2")       U-  =   u   -h  "feu.       or 

el  y 

(2  )      u-3  —   uz  4-  cu. 

These  lead  to  equations  similar  to  (5)  with  the  y -invariants  (arising 
from  (2'))  "being 

h2  -   hz  +■  he  -   d. 


4-   ah   +• 
z 


and  the  z-inyariants  (arising  from  (2  ))  being 


h_  *   c   +■  ch   -  d  . 

3      7  J 


k   ~  e   f  es  -  •, 
3      x  j 


m  ss     c   4-  he   4-  ac   +■  c; 
3      xy     x     7 


It  is  also  apparent  that  whatsTer  results  we  obtain  from  substitution  (2) 
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would  "be  symmetrically  obtained  from  substitutions  (2*)  and  (2»  *), 
hence  at  this  point  We  will  limit  the  majority  of  the  discussion  to  the 
effects  of  substitution  (2)o 

Should  we  be  rewarded  "by  finding  k  =  h,  =  m^  ■  0,  then  (5) 
will  become  a  second  order  hyperbolic  equation  in  only  two  variables 

(7)  u,   +  bu,   +  cu.   +  du-  —  0# 


u,+  b^  +  cu    +  d^ 
jz  z  y 


Ve  may  then  apply  the  Laplace  method  to  (7),  as  described  briefly  in 

18 
Section  I;  or  in  more  detail  in  the  works  of  Darboux. 

Consider  next  the  change  of  variables  (3),  which  leads  to  the 

equation 

(8)      ulz  +  cuf  «  ^  (u)  +  hgux  -I-  ^Uy  +  \xvlz     +     m^  u  y 

where  now  m-^  s  ayz  j-    (ab)2  4-  ca_  +■  abc  «  g.  As  before,  the 

symmetry  of  the  equation  indicates  that  the  following  changes  of  variables 
would  lead  to  equations  similar  to  (8); 

u?  -   V  +  cui , 

J  ■ 

ui  —  u«  +  au^  i 

2  2^  2' 

4   =     \    +  C*2> 

1s\f    Hi  * 

These  substitutions  lead  to  the  respective  invariants 
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5^ 


mj    =      a        +-    (ac)        +■      ha      + 

m      as.     b         4-    (ab)         +      cb„     4-    abc     -    g  . 


m3     ss.    i         4,   (-be)  4,     ab        +  abc 

2  xx  'x  a 


m1      =    c  +■   (ac)  4°    "be        ■+■  abe     - 

3  *y  7  x 


2  ,      1 

m       s    c  4    (Tbc)  4-    ac        +  albc     - 

3  xy  x  y 


2   __ 


In  the  case  of  equation     (8),   if  k    "=■   1      —  hp    =    hl    "=    0,     we  will 
have  reduced  our  system  immediately  to  three  first  order  equations 


ux    +    au     =•     u^ 


(9)  \    +    bul   "     4) 


2     j.    _2    , 


u      +    ci^ 


0 


which  can  he  solved  in  inverse  order  "by  quadratures. 

Thirdly  we  consider  the  change  of  variables  (k) ,   which  leads  to 
the  equation 

(10)        a  ,  4/  au_x  =  ^  (u)   4    l2Uy   f   kguz  4-  m^  u. 

where  m  1  s  d^.  4-  ad  -  g   is  the  first  -x -invariant «,  Once  again 
symmetry  considerations  show  that  the  change  of  variables 

and 

u  -  —  u    4-  au   4  bu   4-  fu. 
-3     sy      y      X      3 
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lead  to  the  -y-invariant  and  -z -Invariant 

m_2  =  ey  +.  eh  -  g} 

m-3  ~  fz  +"  f c  -  g# 
respectively,, 

If  we  are  fortunate  once  more,  such  that  ig  •  ^  =  m^  =■  0, 

we  may  solve  (10)  "by  quadratures,  and  then  apply  the  Laplace  method 

to  the  non-homogeneous  equation  (k) „ 

B,  We  have  in  the  discussion  of  paragraph  A  encountered  some  eighteen 
expressions  which  we  have  labeled  " invar iantsw0  Organizing  these  in 
somewhat  more  orderly  fashion,  we  have  the  following  tahlej 

hg  =  t>z  4-  "be  -  -d 
ho  =  c   4-  "be  -  d 


k  «  a   4-  ac  -  e 
1     ft 

k_  *  c   +  ac  -  e 
3     x 


l-L  =  a_  +•  ab  -  f 
[2  x 


m_  =  a   +  ha_  +  ca  4-  ad  -  g 

j.  yz  z      y  ° 

"2  "  \z   +   a*z  +  Ct)x  *"  *e  -  g 

"3  "  cxy  +  acy  +  tcx  f  cf  "  s 
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■WL  ~   *x  +     ad   "   S 
m-2  s  •_  +•  "be   =   g 

mo   =  f_  4-  cf   -   g 


n^  ="  ayz  4-  (ab)z  4-  ea   4-  abe  -  g 

mf  -  a   -f  (ac)  4*  ba   4-  abe  -  g 

1  yz  y      I 

m_  —  b   4-  (ab)  4-  cb  4-  abe  -  g 

d  XZ                    Z              X 

m3  =  b   4-  (be)  4-  ab  4-  abc  -  g 

2  xx  x    % 

m  »  c   4-  (ac)  *t"  be   4"  abe  -  g 

3  xy  y    x 


m   ^=  c   4-  (be)  f  ac  f 
3     3qr      x    y 


We  now  wish  to  investigate  the  character  of  these  eighteen  "invariants'*, 
as  was  done  in  paragraph  B  of  Section  III  for  systems,  in  order  to  see 
if  the  tform  "invariant"  is  appropriate.  That  is,  we  are  concerned  with 
the  change  in  these  9xpveBB±GD.B   when  we  make  the  various  changes  of 
variables 

(11)  u   —   A  u,  where   A   s  A  (x,y,  z)  is  at  least 

twice  continuously  different iablej 

(12)  x   -  <f  Cx'), 

y  =  <f(y'), 

z  -  Z.U')  3 

(13)  x  » y',    y 
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Let  us  consider  first  the  change  of  variables  (11).  Equation 
(l)  is  then  transformed  into  a  new  equation  of  the  same  type 


(1U)   JfV)  -  «^.,  ♦  a'u^  4-  1,'u^  +  cV   +  i'xi   +  e'i£  + 


hut  with  new  coefficients 


(15) 


1       '    • 
+  f  uz 

4-         '    ' 

T       g   U 

a  + 

1 

X 

it 

■■    h  + 

1 

^x 

c   + 

1 
\ 

}x 
3T 

«  0 


=    d+2.    )x  ^  Ox  ^  1  2  X 

=     a  x   L     ^X    4-  *    ?X     .    1  ^2A 

■b    f  4-   ^      aX     i    a   aX     ,    1  5    X 

X  7x      X  ay      X  5*37 


e  +  ^     -j^r-  ^y    -j^r  jy^-T  ^y-jy  3T^7  ^ray^r  '• 


Using  these  new  coefficients,   let  us  compute  the  three  new     x-invariant, 
k- ,    ]q_,   and  m^   .     We  find 


v.      * 


(16) 


az     4-ac     -e      -    *±  y 

y 


£     =     at     +    aV  -    f '    =     lx 


m!     =    al„  +    h'al    +   c'a^  4-    a  d     -     g    =■  m-i-f . 1  aX  ,  ll   cJA 

t 

The  new  -  x -invariant,  m_-j_>  hecomes 
(17) 


t .  t 


m'  =  d  +  a'd'  -  g'  ~  mol  +  ^  ^X  +  Jk  ^>L  , 

A  7y      \T^ 


p  2*      n- 

while  the  new  mr  invariant,  m^  ,  hecomes 


58 


(18) 


From  gynnaetrj^  similar  results  are  obtained  for  the  other 
h  f   k  ,  1  ,  and  m   invariants  e  Thus  "#e  we©  that  th©  k,  h,   and  JL 
invariants,  which  appear  as  ©©©ffledents  ©f  th©  X,  y,  and  z  derivatives 
of  u,  respectively,  are  indeed  tn»  invariants  under  th®  chang©  of 
variables  (11  J.  The  m  invariants,  ©n  the  other  hand,  do  not  repro- 
duce themselves  ©xaetlj  under  this  ©hang©  of  variables,  but  instead  re- 
produce themselves  plus  a  linear  combination  of  the  true  invariants. 
Hence  we  shall  in  th©  future  refer  to  thee®  m  invariants  as  gmgi_- 
invarlants,  (These  will  appear  again  in  Section  T.) 

If  we  should  make  the  change  of  variables  (12)  and  compute  the 
resulting  invariants,  w©  will  find  that 


■\     «  H^'V 


4>T.V 


(19)  «i  *  qW  V  ial; 

m=l  *  ^x'^y'^z8  «.i. 
Finally  we  not©  that  (13)  merely  interchanges  the  x -invariants 

with  th*  y -invariants^  the  y-invariants  with  the  z -invariants,  and  the 

l -invariants  with  the  x -invariants. 

Ce  Lot  us  now  develop  the  cascade  of  equations  in  th©  manner  of  Laplac® 
and  Darboux,  Consider  first  equation  (5)  in  th©  ©v©nt  that  at  least  on© 
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the  invariants  is  not  zero.     We  wish  to  transform  (5)   into  an  equation 
for    Uj     which  is  in  the  form     (l)   ,   (See  in  this  connection  remark  in 
footnote  (5))       In  order  to  do  so,  we  proceed  to  integrate  equation     (2), 
obtaining 


(20) 


S 


adx 


u    a     e 


\e  u-jtoc    4-    X(y,  z) \j 


where    X(y, z)     is  an  arbitrary  function  of    y     and     z  only.     Then  differ- 
entiating    (20)     with  respect  to    y  and     z  yields 


(21) 


u      =     e 


-J 


adx 


3 


adx 


i.         * 


ux  dx    +  X 


By 


-J 


adx 


ft 


tadx 


ux  dx   +■  X 


(23) 


-s 


adx  ^     Tr    (adx 
(22)         uz=   e      0  ^_  UeJ       u.   dx  *    X 


JF  |1-        -x  -  -    l]     -    £  (J—)  .  "  ^^  ia\dx*j 


Substituting     (20),    (21),   and  (22)   into     (5)     we  obtain 

-Cadx    Tr    Aadx- 
tti      +    bu^    +    cui     4-    dU]_   «■     iie      J  J  uldz        Z 


+  k.   e 


-J 


adx 


-I  adx 
+  1.  e    J 


v*H  -ia 


adx 


(adx 


^dx    +       xj|f 


^adx 


*!  e    °  i|Je         uxdx   +   X^     -      " 


Cadx 


u^dx    4      x|L 


Now     (23)     is  an  equation  for     u-^     alone,  but  it  is  an  integral- 
differential  equation,    interlaced  with  arbitrary  functions  and  their  partjal 


6o 


derivatives,  and  hardly  in  a  manageable  state.  Our  aim  is  to  reduce  this 
to  an  equation  for  u,   in  the  form  of  (l)  and  to  this  end  we  will  dif- 
ferentiate (23)  with  respect  to  x.  Before  we  do  this,  however,  it  will 
he  necessary  to  make  certain  assumptions  regarding  the  invariants  and  the 
coefficients.  If  these  assumptions  are  made,  the  integral  terms  and  the 
arbitrary  functions  will  disappear  when  we  differentiate.  We  will  con- 
sider a  number  of  different  methods. 

METHOD  1 
The  first  assumptions  in  this  method  are 

(a)  l^a  \1    =   n^  i=      0  . 

(b)  a   is  a  function  of  (x)  only. 

Under  these  hypotheses,  (23)  takes  the  form  u^  +■  bu,  +  cis1_m+'  iul  = 


-ladx 


J.^V*  +  x  +  ^{J-Km  4.  ij  *"|[  J.Klfc  4  x] 


\  adx 
Multiply  both  sides  by    JL~,        ,  and  differentiate  with  respect  to  x  to 


get 


ml 


—   6     ^UL  +  *ul  +  cui_>dul)  " 

-  1.  ognl  2_      (Ui      +    bui     +    eun     f  da,)f 
d  x  ax        Tz  -4  -4r         x 

^  adx 

4      \  (UlxyZ  +    tUlXZ  +     CUlxy    +"  *%  *   CXUly  *  Vlz  *    ^KUl)    » 


A  adx  (adx  G 

~       •  ul   4  ul     +     e 


adx 

Ul  ' 
z 
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-  Cadx 
How  we  may  multiply  "both  sides  of  (2k)     "by  m-,©   °     a*id  collect 


terms  to  obtain 


u^    4-  (a  -  Jlog  m,  ) 

X7Z  J7-1    ^Z      "U^Z  +  °%y  +  <%  + 


4  (c(a  -   ^log  n1)  +  c^  -  m^  u^  +  h(a  -  ^log  nj  +  bx  -  m^  +. 

+  d(a  -  ^log  n^)  +  d^.  -  m1)  ux  —   0  . 
Designate  the  coefficients  as  follows: 


al 


a  -   d  log  m-^  b^  m  fc      Cj_  =  c      &,  x  i. 

"J3f 


•l  =    clal   +  C1X     "    *i    j        f  x  =  Dlftl    +  blx     -    m-L  ^ 
«1  "   dlal   4-  «^     -    m-L. 
Then     (25)     is  in  the  same  form  as     (l)i 

^25)       %rz  *"    alUlyz  +   *lU*z  +■  ^lxy  +  dl%+°l^-*-flulx^Slul    **  0.     0 

Our  new    x -invariants  are,   from  the  substitution    u,    s»   Uq     «f-    a  u,, 

^  *    alcl   +■  alz  ~  °i  =  alcl  +   alz  "  alcl     -    %  +  ml  - 

=  alz     "    c0oc  +  ml  ■ 

sc  mi      _    c/  log  m^L     »    cz. 

Thus 

fy     =  k. 


:1  *    (CY  +     ^1°S    mi   \ 
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Also 


\k  -  albl  +  aJ_  -  fl  s  albl-    ^loS  ml  ~  al*i  +  ml   "\  = 
V^x  4-   ^log  "i  )  ^ 

%  -  aldl  +  h.%  +■  C1*L,  +  aiyz   "  %  - 

38  *1*1  "  al^l  "  *1  ^    ^l  "  ^i  #  l°g  *i   "  ci  ^^-°g  "h  . 

_  d    log  m.-^ 

s.  m    -  (  d  *■  b   ^  *log  m1  +  c  J  log  %  •   J  log  \ 

If  now  the  new  x-invariants  are  all  zero,  we  may  proceed  to  reduce 
the  equation  as  described  in  paragraph   A  and  solve.  If,  however,  at 
least  one  of  the  invariants  is  not  zero,  we  are  in  a  dilemma,  for  the 
hypothesis  with  which  we  transformed  (5)  into  (25)  are  no  longer  valid  in 
the  general  ease.  We  can  no  longer  continue  our  chain.  To  enable  us  to  do 
so,  we  must  change  our  original  hypothesis  to  the  following; 

(a)   «!  -  lx  -  "j.  4      °> 


All  coefficients  are  funetiongof  x  only 

(c)    b  ^  c  £L  d  # 
Under  these  hypotheses,  our  invariants  are  also  functions  of  x  only.  W© 
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may  transform  (5)  into  (25)  as  "before,  and  the  coefficients  will  have 
the  same  form.  The  x-invariants  will  he  different,  however,  and  in  fact 
will  he 

k^  =  *!  .  cx   =  *!   -   dI) 
(26)  U  "  ll  -  *x   *  ll   "   dx) 

Hence  k^  =  "]_^  =  m^  ,  our  hypothesis  are  again  satisfied,  and  the  chain 
may  he  continued. 

Observe  that,  as  our  chain  progresses,  our  x-invariants  will  always 
have  the  form 

Thus,  if  our  invariants  are  not  originally  all: xaroj they  can  only  "become 
zero  if  they  are  an  integer  multiple  of  the  partial  derivative  of  d  with 
respect  to  x. 

Note  also,  that  to  have  a  chain  of  y-invariants,  or  z-invariants, 
our  hypotheses  (a),  (h),  and  (c)  must  he  modified  accordingly. 

METHOD  2 


In  this  method  we  assume  first  that 

(a)  kx  =  \     =•  0  m 
Under  this  hypothesis,  equation  (23)  takes  the  form 


6k 


(27) 


Vpadx         r     \  adx 

Ulyz    +  Mlz    +   CUly  +    dul    "    mle  Je  ul 


Multiply  both  sides  "by  © 

ml 
x  to  get 


adx 


,«bc  f    X 


and  differentiate  with  respect  to 


,adx 


e 


(fl,   -       3  log  n^     )  (Uj,       +    bu-^     4-   cu-l     4-    dux)     + 

z  y 


iyg 


X 


£adx 

ml 


h  f     DU., 


xy 


\  +  ^% 


'x"l. 


A  adx 

—       e 


V 


-ladx 
Then  multiplying  by    m    e     v         and  collecting  terms,  we  obtain 

un        +■  (a  -     ^  log  EiNu-i      +  bui       4-  cut       +■  dun      4- 
•'■xyz  i      ■       r J    -^yz  -Lxz  "Tcy  x 

^  x 

(28)  *  jc  (a  -     ^logm^  el    u^    4-  j"b   (a  -       3  log  n^    ^  +  b   1   »i  + 

3  x  J  3  x  J 

4[d(a"   ^LS  )  *  *x   -  ^1  1^*0, 
As  before,  designate  the  coefficients  of  (28)  as  follows g 


a,  «  a  - 


tflog  m. 


bx  -  b^ 


c  ^   d^ ""  cL  ^ 


'1  m  cla!  +  %  ,     f  i  -  Vi  +  blx  ;   81  *  *l"l  +  *lx  "  *L  • 


Then  (28)  is  in  the  same  form  as  (1),  and  is  identical  with  (25) 
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Now,  however,   the  x-invariants  resulting  from  the  substitution 
%*%•+■  alul       axe 

ki^-aici  +  aiz    "    ei  =  aici    -  aici  +  aiz"cix 


■as  a        -     c 
Z  X 


^    log  m-L 


,Z. 


1„     =     a^    ^  aly  -  fx   -  ^  -  hx  -   ^  1|  ^L    } 

=  n^     -     d^.    +    ba    +   ca     +  a^  -  b    <5     log  n^    - 
-  c  d    log  m.3^         _  9    log  n^ 


9  73  x  d$Jd 


x 


If  lfy  »  1^  a  m^  5  0,  ve  may  reduce  as  before,  and  solve. 
If,  however,  at  least  one  of  these  invariants  is  not  zero,  we  are  again 
faced  with  the  dilemma  that  our  original  hypothesis  is  not  satisfied  in 
the  most  general  case,  and  the  chain  cannot  be  continued.  Thus,  we  must 
again  modify  the  hypothesis  in  order  to  iterate  and  continue  the  chain. 
Our  hypotheses  become 

(a)    kx  *  lx  S  0  ; 

(h)    a  =  a(x)  j       b  =■  b(ytz)     c  ^  c^zK 

d  =   d(x)  j       e  ■  e(x,y, z)  =■  ac  j   f  =•  f  (x,y,  z)  »   ab  y 
g  *  g(x). 
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Under  these  hypotheses,  m^ad  -  g,  is  a  function  of  x  only,  as  is 
^log  m-j^     Now  when  we  transform  (l)  into  (25),  the  coefficients 

3  x 


will  "be 
ll 


•a  -   3  log  m1     ,  b-^  —  b  ,  c1  =  c  ,   d^^  *  d, 


el  "  clal  >         fl  -  *lal  >       Si-  dlal  +"dlx  "  %- 
Thus  we  will  find  that 

H  '-     °  »   I4  5  '0  ,   m^  =  m-L  -  dj  , 
hence,  the  hypotheses  (a)  and  (b)  are  "both  satisfied  oace  more,  and 
the  chain  may  "be  continued. 

Observe  that,  as  "before,  as  the  chain  progresses,  the  invariants 
must  have  the  form 

kl+3i    Hll+3i-    °>      ml  +  3i*  ^  "id^  '     i    S  °A,2,...    . 

Thus  if  m^  ^  0,  the  invariants  for  the  Q^  iteration  will  vanish  only 

if  m^  is  an  integer  multiple  of  d^,  I.e.,  m-^j  cL^  for  some  J* 0,1, 2. 

We  also  remark  again,  that  tp  have  a  chain  of  y-invariants  or  a-invariants, 

hypotheses  (a)  and  (b)  must  he  modified  accordingly. 

D.  Now  we  wish  to  consider  methods  for  cascading  the  equations  when  we 

employ  the  substitution  (k),   in  the  event  that  at  least  one  of  the 

•^-invariants  is  not  identically  zero.  For  this  chain  we  wish  to  solve 

(10)  for  u  in  terms  of  u,,^,  and,  after  the  prescribed  differentiations, 

substitute  into  {k)     to  obtain  an  equation  of  the  form  (25)  in  u  i . 

In  order  to  accomplish  this  we  are  again  forced  to  place  certain  conditions 
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on  the  coefficients  and  the  related  invariants.  Thus  we  introduce 

Method  3  s   For  this  we  initially  assume  that 

(a)    12  S   k3  =   0  , 
Solving  (10)  for  u,  we  oh  tain 

u    —     — ■     (u_t     4*  au_-j_) 

Hi        -f  «K- 

and  hence 

by 

H  -  -i    <U-1     +  au-l*  *Ail  )   "  1     9^-l(u-l+au-l)  J 
a-i     ^z  ^2   J7      x 

V-  ^  (U-1^Z  +    aU-ly  +    V-lz  *    az»-ly  *-  V»\l>    " 

~Tf  aair  (u^  +  au^*v*>-|,JJ|  (M^-y+v^ 

Substituting  these  expressions  into  (k)     and  collecting  terms  we  ohtain 
(29)    u  ,    +•  a  ,u  ,  r    +  b  iUj   +•  e  ,u  ,   4  i_iu  _  +  e  u  -,  4. 

"*  f-lu-lz  +  s-lu-l  *  °J 
where   a-X  =  a  ,   h_x  =  "b  -   <y  log  m-l  ,  c^  «*c  -  J  log  m-l 


7 


68 


d-_l  —  &  -  c  a  log  m^  -  b  <$  log  m^  -   a^Lcg  m_i 

3  y       ~5~£~~     "5  j^z 

+    °g  m^   ^  i0g  a  1 

TV"'"      "7T^  ' 

6-l*   a-lZ    +    a-lC-l   *      f"l   ^-ly   +a-l^-l   >      aJld 

8"lr:  a"1yz  +b-la-lz4c-la-l5r',-a^l  "  m-l  ' 

and  (29)  is  of  the  desired  form.  Computing  the  new  -x-invariants,  we 
obtain 


k.   =r   C    —   a   " 

6     x      z 


^:x:^y 

^iog  m-l  - 

"  ■  — ■■  ■  —  ■  •  * 


and 


m  .  -=  m^  -f  d^  -  O       (c  dlog  m_x  ^  h  J  log  m^   ^*log  mml  „ 


T* 


7~7 


_   d  log  m^   ^logm-]L)  —  ay  (c  -  ^logm.j) 
~   az  (h  -   3  log  m^)  +  a^  # 

c)y 

As  in  paragraph  C,  we  will  he  in  a  dilemma  if  the  new  -  x- 
invariants  are  not  all  zero,  fcsrthen  our  initial  hypothesis  (a)  will  not 
he  satisfied  in  general,  "by  (29),  and  our  chain  cannot  continue.  To 
escape  this  trap,  we  once  more  impose  additional  restrictions  on  our  coef- 
ficients. Our  revised  hypotheses  will  he 
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(a)  12  »  k3     S    0    • 

(b)  a  =  a(x),         b  ^b(y,z)j        c   «c(y,z)  , 

d    =  a(x),       ©    =-6(2,7,2)    =  a(x)c(y,z)   , 
f  -f(x,y,z)    =  a(x)b(y,z)  ,     g  ag(x)  . 

Under  hypotheses     (a)     and     (b),     m_j_    "becomes  a  function  of    x 
only,  and  hence  the  coefficients  of     (29)     "become 

a-lm  aj         k.i^  b,       c„-|_   —  c,       d  a  =    cL>       o„i  ~    e?     '•!  —  f 

8-1  —  ad  -  m_-j_  , 
With  these  coefficients,  the  -x -invariants  for  (29)  are 

I5  S  k6  *  °>    m-4  *  m-l  +"  *x  > 
which  means  our  hypotheses  (a)  and  (b)  are  again  satisfied,  and  the 
chain  may  be  continued,,  As  observed  in  methods  1,  and  2,  the  invar- 
iant m.-jL  _  o*  can  vanish  only  if  m„i  is  an  integer  multiple  of  d^, 
that  is,  m_-j_  ■  -Jdx.  We  may  add  at  this  point,  that  if  m  is  a  positive 
multiple  of  d^.,  we  should  use  substitution  (2),  while  if  it  is  a  nega- 
tive multiple  of  d^.,  substitution  (k)     would  be  the  more  advantageous. 
One  more  it  is  worthwhile  to  note  that  if  we  desire  a  chain  of  -  y  - 
invariants,  or  -z-  invariants,  hypotheses  (a)  and  (b)  must  be  modified 
accordingly c 

E.  Finally  we  consider  a  method  for  cascading  the  equations  when  we  employ 
the  substitution  (3),   in  the  event  that  at  least  one  of  the  invariants 
in  (8)  is  not  zero*  Although  this  substitution  appears  the  most  natural 
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in  that  our  equation  (l)  is  immediately  "broken  down  into  three  first 
order  equations,  we  will  find  that  the  most  stringent  conditions  on  the 
coefficients  are  required  in  this  case  in  order  to  generate  the  chain. 

Method  k:       Our  initial  assumption  in  this  method  is  that 
(a)  hp—  k,  SI,  s  0,  We  must  consider  the  system  of  equations 

u^  +  au  "»  u^ 

\   +  tul  =  Ul  J 

2   .    2  _   2 
ulz  +  cul  ~  mlu  ♦ 

If  we  solve  the  second  equation  for  u^,  substitute  this  value 
into  the  first  equation  and  solve  for  u,  and  then  substitute  the  result- 
ing expression  into  the  third  equation,  we  obtain 

2       2     2   -iadxf,  Jadx-bdy  r  jhdy  2 
ux   4-  cui  =  mi©         \  e        1  e     u  dy  +" 


4yU,z)      to  +  x(y,z) 


where     Y    and    X     are  arbitrary  functions  of  their  respective  arguments. 
Ve  then  proceed  as  "before: 

W 


e 
m 


fadx      T    p  2\       r  .\adx-bdy  f  r    Jjbdy  2  7 

_    R  +  cui  \sif  H9       ^  +  TJte4-X. 


1 
Taking  the  partial  derivative  of  "both  sides  with  respect  to    x, 

(30)        _e__       (a  -      £log  m-L  )(U1     +  cy?   )  +  u2      +     2         c  u?  U 

4     i         s  x        z      x      lxz     **     x  J 


Jadx-bdyf  f  Jbdy  ) 

I  Je        ul>^    *  YJ 


-     e 
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Ybiy-adx 
We  then  multiply  both  Bides  ©f  (30)  by     ©  u  and  take  the  partial 

derivative  with  respect  to    y,     to  obtain 

(31)  n2p22  \  2       2  2  2. 

+  (■7     -    d    logiLjKttJ   fCW1)4   (a  -      ^loffl^Huj^    +  «i->CjH1) + 


d  *3  7  T 


jh&y 


2  2  2  2  2l  AW     2 

V+  %  +  VS.  +  '*%•  +  «ur*i]  *   * d       ui . 

0    \-bdy 

Finally,  multiplying  (31)     by    mj©^        ,  and  collecting  terms, 


we  obtain  an  equation  for    Ui[,     which  is  of  th®  same  form  as     (l),  hut 
whose  coefficients  are     a_  2  a  -      ^iogm,    j     b£  *  b  -    plogm-j     c1  =   c  « 

^2         2  2.2  222,2  -2  _  a2k2    .  a2      . 

dl  *  alCl  4  Cly*       ®1  *  alel  +  %   ^    *1   *  Tl    *     ly    ' 

2         2,2  2       2. 2  2  2  2  2  2 

^s  y^fapl^  4-  (aiei)y4»c.         -    mx  , 

Computing  the  new  invariants,  w©  find 

2  2  32    _   .2 


lz  ^  »!•!  -     d!   -  -hlz  -  ely  ^ 


*;_  + 


2       22       2       2       2 

H   s    *1.  +  alcl  "     @1   *  alz  "  clx:> 


2 


A.?-  f? 


2  £  ,    ,  2. 2,      ,        2  2.     2,  2  2  2 

"\  =      alyz      ■•"  (albl'a+       Claly   ■**  altlel       ■      aS 

?  2  2_  2  2^2 

-  »i  +  av  +  wi>,  +aiyci  4°l3fy  • 

We  are  faced  with  th©  same  difficulty  as  "before.  In  general,  the 
invariants  h«  and  lq.  ^-H  not   vanish,  henc©  it  will  b©  impossible  to 
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continue  the  chain  in  this  manner „  Wo  will  "be  able  to  continue  the  chain 
if  we  revise  our  hypothesis  to  he  (a)   hg  S-  k,  =1.=  0 

(h)   All  coefficients  are  functions  of 
j     only,  except  c,  and  that  c 
is  a  constant. 
With  these  hypotheses  our  new  coefficients  "become 

al  =  a|  *1   *  '   ^lo&^>     ci-  c  °>     di'^lcl  °> 

c)  y 

ei  =  44°>  A 9  ■M*  ai7>  &  ■  ^M  -  ^  >  ** the 

new  invariants  became 

*5  s  **  s  h  s  °j 

m?  -  ml   ■+■   Cay  • 

Thus  the  hypotheses  (a)  and  (h)  are  satisfied  once  more,  and  the  chain 

may  he  continuedo  We  observe  that  in  this  method  the  chain  will  terminate 

2 

with  all  invariants  zero  if  and  only  if  nu   is  a  negative  integer  multiple 

of  cay. 

These  methods  described  in  paragraphs  C,  D,  and  E  are  not, of 

course,  the  only  methods  by  which  chains  of  equations  could  be  generated. 

We  could  do  a  method  similar  to  method  1,  say,  with  k_  3  0  and  J,  s  nu? 

or  we  could  do  a  method  similar  to  method  2  with  k^  the  non-zero  invariant, 
or  lj_  the  non-zero  invariant.  There  are  many  variations  on  the  theme 
presented  here,  but  it  is  the  theme  itself  which  is  the  important  thing. 
In  any  such  method  certain  restrictions  of  the  coefficients  are  necessary 
to  iterate  once,  and  other  restrictions  are  necessary  to  iterate  more  than 
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emcee  The  "basic  idea  is  to  find  an  expression  for  u  in  terms  of  the  new 
dependent  variable,  which  may  or  may  not  he  an  integral  expression*  We 
then  differentiate  this  expression  as  required  hy  the  particular  equation, 
i.e.  the  equation  involving  the  invariant  coefficients.  Substituting  the 
resulting  expressions  into  this  equation,  we  reduce  it  to  an  equation  in 
the  new  dependant  variable  alone;  if  necessary,  we  separate  these  ex- 
pressions from  the  invariant  coefficients  and  differentiate  again  suf- 
ficiently to  remove  any  Integral  signs  which  may  appear.  If  we  can 
accomplish  this,  we  can  in  general  iterate  the  equation. 

F.  We  consider  now  the  most  general  form  of  (l)  which  can  he  reduced  hy 

these  methods  described  in  paragraphs  C,  D,  and  E.  The  question  is,  what 

must  the  coefficients  of  (l)  he,  in  order  that  (l)  may  he  reduced  to  a 

second  order  equation  hy  these  methods,  when  the  x-invariants  or  the  *x- 

invariants  are  not  all  originally  zero?  We  observe  first,  in  methods  1,  2, 

and  3>  that  in  order  for  the  m-invariants  to  vanish  after  n-1  iterations 

we  must  have 

g  -  ad  =  nd^  ,  where  n  =■  ±1,  ±2,  ... 

or 

(32)  a  4  ad  -  g_ 

n     n 

Equation  (32)  is  integrable,  and  integration  gives  us 


■S     MeS 


dx 
(33)  dseO  le  ^  g  dx      plua  an  arbitrary  constant 

U  n 

which  we  choose  to  be  zero. 

Using  (33)  and  the  hypotheses  of  method  1,  we  see  that  our  original 
equation  must  be 


7^ 


(3*) 


V+a(x}v4  0        r      s&(^%+Mf 


-4-  8(x)(1»y+\+  «)  ■  0    ,   nsl,2,  0o 


for  method  1  to  "be  applicable 0  Using  (33)  and  the  hypotheses  of  methods 
2  and  3j>  we  see  that  ©up  original  equation  must  "be 

(35) 

_j.  a(x)e(y,  z)"UL  +  a(x)b(y,  s)«z  f  g(x)"^  s  0  ,   n  *  1,2,  <,<,<» 

for  methods  2  or  3  to  be  applicable. 

Finally  we  consider  the  form  of  the  coefficients  when  method  k  will 
be  applicable o     In  this  case  we  observe  that  in  order  for  the    m    invariant 
to  vanish  after  n-1  iterations,  we  must  have 


g  =  abc   ss.  nca_     ,     n  «•  1,  2, 

or 

(36)  a      +    ha    -    g_   # 

*  n  nc 


o  o  o 


Equation  (36)  may  be  integrated  to  give 
(37)  a  =     e       On      \  ©  Jn      gjj  # 


c  n 


Using  (37)  and  the  hypotheses  of  method  \7   we  see  that  our  original  equation, 
must  he 


«       +  («      Jn       \e   Jn    ^^yju      +  h(y)u     4- cu       4- 

-  Chdy  r    (bdy  * 

+  Cb(y)ux  +  (ce      Jn     \®Jn      _£ay\iUy  + 

(38)  «J  en/ 


c   : 
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+  /JL   +  MnIi)Q    -  C^  fe  i^JL   dy  )  \  +  s(y)u  =  o j 


for  Method  lj-  to  "be  applicable. 


SFOTIOH  Y 

i 

HYPERBOLIC  EQUATIONS  OF  N™  ORDER 
IN  I  INDEPENDENT  VARIABLES 

Ao  Having  considered  the  extension  of  the  Laplace  cascade  method  to  the 
equation  of  third  order  in  three  variables  of  the  mixed  derivatives  type, 
we  will  now  turn  ©ur  attention  to  a  generalisation  ©f  these  results  •  That 
is,  we  consider  an  equation  of  n*k  order  in  n  independent  variables, 
in  which  only  the  mixed  derivatives  appear,  and  we  shall  attempt  to  show 
how  these  methods  can  he  applied  to  such  an  equation.. 

Before  proceeding  to  the  n**1  order  equation  we  will  first  intro- 
duce an  operator  notation  which  will  simplify  somewhat  the  invariant  ex- 
pressions and  relations,  which  tend  to  become  awkward  as  we  increase  n. 

In  order  to  do  so,  we  mention  here  a  rule  which  has  been  devised  by  J.  Le 

19 
Roux   and  shall  hereafter  be  referred  to  as  Le  Roux's  rules 

"For  any  differential  operator  D  which  is  a  polynomial  in 

and  Q      the  following  relation  holds  s 


D(uv)  ^uP(v)  f  9u  D^v)  •  §u  Dy(v)^2j 
dx       £y 


22u  D12,(v)+2^2u  dI1(v)j-  ^Ud"  (vJ  + 


-] 


where  D°  s  *D   ,  ,  etc, 
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This  rule  is  an  extension  of  Taylor's  formula,  and  is,  of  course, 
completely  analogous  to  Leibnitz'  Rule  for  the  ntJl  partial  derivative  of 
a  product  u  v.  It  can  he  similarly  extended  to  products  of  more  than  two 
functions.  As  an  example  of  this  rule,  consider  the  operator  D(u)  — 
=  u +-  au  +■  biL,  +  cu  —  0.  Here  D  =■  ^  d     4_a^  j.  b  ^  -HC,  hence 

D=d  +  a  ,  I)'  »  d  H  |  D"  =  I,  the  identity  operator. 

<5F  jx 

For  this  operator  D,  it  is  a  simple  matter  to  verify  the  identities 

(1)  DV(u)  a  D(u)  -h  hu 

y  x     ,       ; 

(2)  D^u)  =  D(u)  +  ku^ 

where  hsa^-t-ab-c,  k  »  b  4-ah-c   are  the  two  Darboux  invariants, 
^x  y 

Then  denoting  u-^sD^u)  and  observing  that  D(u)  =-  0,  we  have  from  (l) 

(3)  Dy^)  *  hu  , 

Operating  on  (3)  with  J>x,   we  obtain 

w      Vy(ux)  =-D;(hu). 

Using  the  identity  (2)  to  evaluate  the  left-hand  side  of  (k)   and  Le  Hoax's 
rule  to  evaluate  the  right-hand  side,  we  find 


or 


sz    h  Ut  4-  0  h  u 

57 

=■  h  Un  4  Jh  .  h  .u 
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D(u-l)  -   Jlogh  Dy  (vlj)    +  (k  -  h)ux  -  0  . 

This  may  "be  written  in  the  form  D,(u,)  "=  0 
vhere    Dn  -   ^  d     x  a,  c        +   hn  »   j_  c_  , 

^77    J*     dT 


a^—  a  - 


h-»3 


^logh 
47 


i±  =  c  -  h  )logh  ^.hy  -  ax, 


20 


77 
and  these  coefficients  agree  exactly  with  those  ohtained  hy  Darhoux' 
for  the  first  equation  of  the  cascade  method. 

Our  prohlem  now  is  to  expand  these  procedures  to  the  n*11  order 
equation  in  n  independent  variables,  of  the  mixed  derivative  type,.  Can 
we  in  this  manner  predict  the  form  of  the  "invariants",  the  corresponding 
identities,  and  methods  for  cascading  the  equations? 

B.  Let  us,  for  convenience  sake,  change  the  notation  slightly.  Instead 

of  indicating  the  independent  variables  hy  x,  y,  z,  ...„,  we  will  find  it 

th 
more  advantageous  to  label  them  x^,  Xg, ...  .x^e  Then  let  the  n   order 

operator  with  which  we  will  be  concerned  he  defined  hy 


(5) 


D(u).   >_  au  ,,!  ui  i  ..  .i  ,  i*  =  o,  1,  2,  „<,•., n 


where  we  always  require  that  i*  4=  1.     ^     and  that  ij-  ij^-l   *%  aa<i  o&ly  if 
^J  +  l*00     T^ae  BVm  is  *°  te  'talc®31  OFer  a11  acceptable  combinations  of  the 
ij  where  i^^o,,,  •<,,&<>     In  (?)  v©  vrite    uj«  jjfy  >  <J~°>  **   •••»*  a3Qd  &*fijas> 
*u    S  Uo         Finally  we  require  that  in  all  of  these  operators  we  will  set 


xo 
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a.       ••»£  —  !•   To  illustrate,  let  a  =■  2.  Then 
For  n=  3  wn 


+  a001  i!L  +  a002  ^L  -I-  a003  ^L  "^aOQQu  f 

and  the  brevity  of  (5)  "becomes  more  apparent  with  increasing  n.  Having 
seen  these  examples  in  explanation,  there  will  he  no  loss  of  clarity  if 
we  drop  the  zero  sub  scripts  of  the  a's  wherever  they  may  appear.  Thus 
for  n  =.2,  we  write  D(u)  =  u_p  4-  ^gUg  +  a-^  -f  **• 

In  order  to  point  the  way,  let  us  briefly  summarize  the  2nd-  order 
results  of  X&rboux  and  our  3   order  results  in  the  light  of  this  notation* 
JBgr  the  second  order  equation  we  have  the  following  operators: 

*>  -    L    *      +  al  *   +  a2L  +   a 

axx  pg        j£[        ^  > 

1    ax-       * 


2 


i 


^  =  i 


The  two  Darboux  invariants  are  given  by 
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b  —  0     a.,  4-  a  a.  -  a  =  D  (a_)  -  a 
t— r  1^21       2  1     > 

k  -   5  a2  4  a2a2-  a  =  D^)  -  a } 
0*2 

while  th©  identities  (l)  and  (2)  became,  respectively 

D^d'Au)   "■  D(u)4   hu  =.  D(u)  4-  \  T>2(&i)   -  al  u 

D-Jp^Cu)  ~  D(u)  +.  ku  ^  D(u)  +  (^(ag)  -  ai  u. 
For  the  3   order  equation,  the  operators  are  as  follows: 


D 


<L<L    d    fa^      i    il+a^^^.  +023    5      d      + 

0X1  d*2    dx3  ^*1^2  d*l  5*3  Jxg   <5)x3 


+  a^  +  a^2  +  a3^  +  % 


D^  =  ^     d 


44+ai2<^ai3fe*aij 

D2  s  »$—   +"  al2    ^       +•  a23    ^       -*-   a2 

D3=LL    +  a  ™    *      +   a23  i_  +   a3  s 

D12  =  j~-  4-  al2  5 

D*'   =    d      +  a 
13      S3£  13-> 

D23  "  -2L  4-  *23  * 

and  associated  with  this  equation  are  the  invariants 
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h3  =  B^U^)-^  4  -  D^(a23)-a 


u 


k3  ~  B&(a12)-aS  "I  "  B23I)]2(al3)-a 

Xl  =  ^JS*-^  ™3  =  ^i'V"8- 

1  -=.   d"  (a  )-a  m2  =  D"D"(a  )-a 

2     23N  13   3  3     23  13  12; 

"l  =  Dl  (a23)-a         -        ^l"   ^^l^ 

"fe  -  D£  (al3)-a  m-2=   D^)"a 

n^  =-  D3  (a12)-a  m-3  =   Dj^8^)"* 

A  typical  Identity  involving  these  invariants  is 

t  it  "  H  * 

DjD^U)  -  D(u)  4-  (D12(a23)-a2)u2+(D13(a23)-a3)u3-^(D1(a23)-a)u# 


C,  Now  let  us  generalize  these  results  to  the  n**1  order  linear  differ- 
ential equation  of  the  form  (5)o  In  order  to  demonstrate  the  magnitude  of 
the  problem,  let  us  first  compute  W(n),  the  number  of  different  identities 
associated  with  the  n  ^  order  equation,  D(u)  =  0. 

An  identity  for  D(u)  is  a  combination  of  differentiated  operators 
D*,  D  ,oo.  such  that  successive  application  of  these  primed  operators  to  a 
function  u  =  u(x),  where  J-Cx^Zg,  ...,3^),  will  yield  D(u)  plus  a  linear 
combination  of  the  mixed  derivatives  of  u  of  order  less  than  n.  We  shall 
denote  such  a  linear  combination  "by  R,  for  remainder,  and  define  the  coef- 
ficients of  R  to  he  the  invariant  coefficients  (or  simply  invariants) 
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associated  with  the  operator  D.  We  shall  show  presently  that  the  order 

of  B  must  he  less  than  n-1,  although  this  is  not  evident  a  priori. 

To  he  a  systematic  we  proceed  as  follows: 

i 
Select  arbitrarily  a  primed  operator  having  order  n-1,  i.e.  D. 

for  same  i«l,2,  ...,n.  There  are,  of  course,  n  choices  for  this  oper- 
ator. To  obtain  an  identity  as  defined  above,  there  is  only  one  choice 

of  operator  D    such  that  D  D*$)F©(u)+R.  This  operator  must  he 

(n-1)  ^ 

Dn  0   .  1    .     .         ,   since  this  contains  the  differentiation  O     which 

t 

is  missing  from  D..  Therefore  we  have  exactly  n  identities  of  the  form 

t  (n-l)  t    x 

D±  Di  ...  i-1  ifl...nW  "  !>(«)  +  H, 

i^^  jLfdy  . .  o,n. 
Next  select  arbitrarily  a  primed  operator  of  order  n-2,  D   fa** 

some  i  =  l,2,  ...n  and  some  J»l,2,  ...,n,  J^i.  There  are  n  choices 

"    i 

for  i,  and  n-1  choices  for  J,  hut  since  D  *  D  ,  we  total  only 

ij       Ji 

&(a«l)  -  |  p\  such  operators,  where  (7\  denotes  the  appropriate  "binomial 
coefficient.     If  we  wish  to  obtain  an  identity  when  we  apply  such  an  oper- 
ator, we  must  apply  it  to  one  of  three things: 

1)  DU"2)  <»> 
x>       vl  2  ...  i-1  141  ...  J-l   j+1  ...n       . 

2)  D(n-1}  D(n-1}  (u) 

12   ...i-1  i4l  ...n  1  2   ...  J-l    j+1  ...n 

o)     Jn'V  D.(n-1)  (U) 

5}     vl  2   ...J-l  J41  ...n  1  2   ...i-rl  i41  ...n 

Therefore  we  have     3H  identities  of  the  form 


8* 
D   L(u)  b  D(u)  +  R 

whore  L(u)  is  on©  of  th©  three  expression  shown  above. 

If  we  continue  to  count  th©  identities  in  this  manner,  we  find 

that 

n-1 
(6)         N(n)  »   ]T  &J  (*)  , 

where 

0lj  "l^ivpi*  t*1®  B**  of  aH  ordered;  integral  partitions  of  J, 


^"fv"*^:*!!^] 


iS.^™-L,  ^,  000  j      J 


To  verify  this  formula  it  is  only  necesary  to  count,  as  we  have 
already  done  for  the  terms  J^l  and  j=2,  the  identities  for  the  j 
term.  For  each  of  these  identities,  the  last  operator  to  "be  applied  is 
D.  .      .  ,  where  the  i   form  all  the  possible  combinations  of  n 

1*  ip     Oil     ll  P 

integers  taken  j  at  a  time,,  Hence  there  are  (n)  choice  for  this  oper- 
ator. To  complete  an  identity  we  must  apply  this  operator  to  a  combina- 
tion of  operators  which  contains  each  of  the  differentiations  & 


p  s  1,2,  ...,  j,  once  and  only  once.  Thus,  for  example,  we  may  have  a 
combination  of  j  operators,  each  of  which  has  on©  differentiation  O       \ 

or  we  may  have  a  combination  of  J-l  operators,  on»  of  which  contains 
th©  sefcood  order  differentiation  O \_  ii£  while  th©  others  ectafcalni' 


only  first  order  differentiations  d_     j  and  so  forth,,  The  nuoiber  of 


1P 


choices  for  each  operator  can  only  "be  the  number  of  combinations  of  the 
differentiations  still  to  be  included,  taken  as  many  at  a  time  as  the 
order  of  the  operator  "being  considered.  Finally,  since  the  order  in 
which  these  operators  are  displayed  is  a  factor  in  determining  which 
identity  we  have,  we  must  consider  an  ordered  integral  partitions  of  J. 
The  summation  in  (6)  is  taken  only  up  to  n-1  since  for  <J  —n,  we  have 
the  trivial  identity  D(u)= D(u),  and  this  is  not  to  Tbe  counted  in  our 
determinations.  If  we  tabulate  the  appropriate  values  for  n=2,3,*i-,  an^L 
5,  we  have 

a=   1      1(2) »  2 

a2=  3      N(3)=12 

a  a   13      *W  -  7^ 
3 

a^  s  75     1(5)  =  5^0 

D.  We  turn  now  to  a  consideration  of  the  invariant  character  of  the 
coefficients  of  an  arbitrary  R.  Consider  first  the  identities  of  the 
form 

"ftlT..i*i+i...»,,iw~JW  +  B- 

The  expression  d'  (u)  is  of  order  n-1  and,  in  fact, 


y   tfi)  it 


8.5 
n  n  ^  v 

'  "o"  ^-  A-  0.   ;.iui...ln>.i   »v. ^_     U  + 

a** 


We  may  writ©  this  more  compactly  as 

(8)  ^H^t^fcMS  ±  w&2j*---*i«. 

To  illuminate  the  notation,  the  superscript  on  u  Indicates  the  order  of 
differentiation,  the  subscript  on  u,  which  is  in  parentheses,  indicates 
the  variables  with  respect  to  which  u  is  not  differentiated  and  the 
subscripts  of  the  a  coefficients,  which  are  in  parentheses,  indicate 
the  integer  or  integers  omitted  from  the  subscripts  of  the  a  coefficients 
in  D,  For  example, 

a/^\=L  ai2  oo.  j-1  ,J+l...n 

a/.%  \""~a   ooc  j -»l  j<|-1. . o  1c— 1  xc4*l»e»n* 

Observe  that  a,   ,  =  a,  )f   hence  the  factor  h     far  similar  reasons,  the 
q    term  in  this  expression  contains  the  factor    1   « 

r^in 

If  we  scrutinize  equation  (8)  we  find  that  there  is  just  one  term 
of  order  n-1,  and  its  coefficient  is  1;  there  are  n-1  terms  of  order 
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n-2,  and  their  eoef  fie  lent  s  are  the  a/*\*  there  are  C31^1)  terms  of 
order  (n-3),  with  coefficients  a/jk\j  and  .in  general  there  are  (JJ-a-0.) 
terms  of  order  n-q..  Consider  now  the  effect  of  the  operator 
^...i-l  i+l  ...  n  ~   +  a(i)>  ^^  tt  is  applied  to  DjU).. 


n  /_  ,\       n 


»<**)  »;  (U) .  u<n> «.  y  a,  .■J"f)  tI[k,u         \  ,n-2)  + 

i...i.i  ua..#  iv  jil   (J)  (J)    +jti  I  JjJ*"fc)"fr)hu>  + 

3*1       k*i  .  dXi 

J*k' 

Equation  (9)  tells  us  immediately  that  there  are  no  invariant  coefficients 
for  the  terms  of  order  n-1.  This  statement  is  readily  seen  to  he  true 
no  matter  which  identity  we  consider.  Thus  a  discussion  of  the  invariants 
must  "begin  with  the  terms  of  order  n-2. 

From  the  fact  that  the  indices  1,   J,  and  k  are  distinct  we  con- 
clude that  W  ij )  ^  u(  jk)  f01"  ^  «5»  We  further  note  that  (9)  contains 

(n-l)  +•     ( n-1 ) (n-2)  ~(o]   different  mixed  derivatives  of  order  n-2,  and  ': 
2      w 

these  are  all  the  possible  mixed  derivatives  of  this  order.  Thus  in  order 

to  ohtain  the  identity 

(n-l)  *  /  \    /  \  i  « 

Dl    i-1  1+1     .»-<*)=»(»)  +  *> 

i.e.  to  make  the  expression  D(u)  appear  on  the  right  hand  side  of  (9)>  we 

must  (at  the  very  least)  add  and  substract  to  the  right  hand  side  of  (9) 

the  terms  a/  **5  n»   ^als  results  in  the  invariant  coefficient  of  order 

n«£: 


8$ 


(10) 


N. 


(jjHji  +  au)  in  "a(jD  '  w3*2*—*  **J- 


Sine©  there  are  n  choices  for  i,  and  n-1  choices  for  j,  we  have 
a  total  of  n((n-l)  invariant  coefficients  for  the  terms  of  order  n-2. 
We  will  now  prove  the  following: 

Theorem  7:  The  invariant  coefficients  for  the  terms  v-i**)   of  order 

n-2  which  appear  in  the  remainder  R  of  any  identity  are  of  the  form 


^a(j)  +  aQ)a(i)  "a(ji) 

3^T 


These  "invariants  are  all  true  invariants,   i.e.  invariant  under  the 

change  of  variables     tt  =  Au,   and  no  invariant  coefficient  for  terms  of 

if  tt 

order  less  than    n-2     ean^in  general  he  a^.-true  invariant. 

Consider  any  identity  given  "by  a  sequence  of  operators, 


.(a) 


S» 


•  B 


W )    >> 


(u)  =D(u)  +■  R. 


vit.0o±&     J|»»«Jfc    k-L*..*^  11»*--1b 
where  (n-a)-f-  (n-Ta)  +  ...  +■  (n-s)  ss  n, 
If  we  write  this  out  in  more  detail,  we  find  that  such  an  identity  is 


actually 


n 


^(n-a-Jj 


♦  +£ 


(h-a-a) 


»*•      5xi  •  •  •  JX^       jml      »J£  .  .^Xl . .  .^xn         ^  iK=l  ^x  ...^..^ 


,(n-b) 


yn-b-1) 


T  ^^n"s)  „  #-  n(h-s-1) 


*  •  • 


=  D(u)  4-  R  . 
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By  an  inspection  of  the  left  hand  aid©  of  (ll)  we  can  determine  how  each 
term  in  D(u)  is  mad©  to  app©ar  on  th®  right  hand  side  of  (ll).  To 
oTbtain  th©  n   order  term    a   u  s   Wi  must  combine  th®  highest  order 

differentiation  in  each  operator  D>*'    ,  D^  .  /^i    u  >   smd  ^P1^ 

,  _*        *1##***   Jl"-JV  Ji£i».»K 

th®  result  to     ^     U      .  ^he  determination  of  this  nth   order 

0  3C"i  o  o  o  Gp^Ci  •  •  a  Cf  %n 

(i^  ll««o,l  ) 
J"       B 

term  is,  in  fact,  th©  criteria  used  to  determine  which  combination  of 

operators  yields  an  identity.)  The  terms  of  order  n-1  which  arise  are 

obtained  in  two  ways0  They  are  obtained  by  combining  th®  highest  order 

differentiation  in  each  operator  D^;,,!  ,  J>\   J -4j  y^i  ,6ok  9   a*14  aPP^= 

x  a    b    j_    q^ 

ing  the  result  to  ©aeh  term  of  the  simi 


(12) 

(i,j  »l1«..l83  i*  d) 

Using  Leibnitz's  rule  w®  see  that,   among  others,  the  terms 

_  (n-l)  o       .   0 

a(j)    *  (j)  «3  =  l,2,oo. ,n 

must  appear o  Such  terms  of  order  (n-l)  also  arise  when  we  combine,  on© 

by  on®,  a  differentiation  of  second  highest  order  in  any  one  operator 

•*    ,  ,    dI^  >  with  th©  highest  order  differentiation  in  every  other 

i.  oooi    ...J  K-^c.eJtq  ■»- 


operator  ^U^mi^Jt^i 
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5  (n"s)u 

i  #  lx*««lfl 

The  coefficients  of  the  terms  of  order  n-2,  -which  are  the  terms 
we  are  concerned,  with  in  this  theorem,  arise  in  four  ways.  We  may  ob- 
tain such  terms  "by  combining  the  highest  order  differentiation  in  each 

operator  D^    ^   ...  D^  . .  .k  and  applying  this  result  to  each  term 

x        a  a,         *i 

of  the  sum  (12).     Again  using  Leibnitz's  rule,  we  see  that,  among  others, 
the  terms 

9  a(j)     .   u       f  i,J  =  1,2,. ..nj   i  ¥  J; 

-J—  (ij) 

1»JT  l-i  o  •  .1 

must  appear.  If  we  combine,  one  by  one  a  differentiation  of  second 
highest  order  in  any  one  operator  D*a'       D*4'  .  ,  with  the  highest 

order  differentiation  in  every  other  operator  D*a'  4-w..,d£5)  .   and 

i2»o.ia>   «»k^...k 

aPPly  this  result  to  each  term  of  the  sum  (12),  we  see  that,  among  others, 

the  terms 

(n-2) 
a(j)  a(i)  U(U)         *»J  *l»2,...,nj  i  +  J 

i,j  +  V#1a  ) 
will  arise.     In  a  similar  manner,   if  we  combine  a  differentiation  of  third 
highest  order  in  any  one  operator    D^  0,o±  ,  •••Dk!r...k      with  the  highest 
order  differentiation  in  every  other  operator  Dj^i  r°jc  ...k      and  ^P1^ 

this  result  to  (13),  we  obtain  the  terms  of  the  form 

(n-2) 

a(Jk)  U(jk)     ' 

while  if  we  combine  the  highest  order  differentiation  of  each  operator 


90 


d!         «  &k:        v  >     smd  ^P1^  tJlis  3P»snlt  to  each  term  of 


(n-s-2) 

z 

<3*k       (±,J.>k»l....l  |  i*j,k) 


a/,i.\    *       U 


di^TTryx^TTry^ 


we  must  again  obtain  terms  of  the  form 

(n-2) 

Comparing  the  coefficients  of  these  terms  of  order  (n-2)  and  recalling 

the  procedure  necessary  to  make  D(u)  appear  on  the  right  hand  side  of 

(9),   we  see  Immediately  that  the  invariant  coefficients  for  the  terms 

a)*I%  *  which  appear  in  the  remainder  E  must  he  of  the  prescribed  form 

(10). 

Knowing  that  these  invariants  must  he  of  this  form,  we  proceed  to 

show  that  they  are  true  invariants.  We  consider  the  equation  D(u)-0 

and  the  change  of  variables  u  «  Xu  f   where  X*-  X(5)  is  not  zero.  Using 

the  chain  rule  for  differentiation  of  a  product,  and  employing  the  notation 

of  (8)  we  have 

n  n      n 

u(n)^Xu(nM^X  $**\   ill   ?x     z<*f+ 


(1*) 


(n-1     \_(n-l)      JL    nv      _(n-2)  *l    Ji     \2  \       4».3)    , 

W             v  ;      j=a  3Z,    lw  2    M  k«l  a  1 ■  d xte    (ijk) 

J*i        J  J#i  k*l        *     '* 

<i*k  i*l  P  n 

(n-2)     v  _(n«-2     3    n*        U-3).  ,    ._,   0           _     ..  , 

k^ 
k*i 
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and  similar  other  expressions  involving  derivatives  of  order  lass  than 
(n-2)o  Substituting  the  relations  (ik)   into  the  equation  D(u)  »  0,  and 
dividing  through  "by  A  ,  we  obtain 

5(5)  =    uU)  +  f   (a,.  +   )log\  )  u(;;f  + 

i_     (i)    -«; —     (i) 


1*1 


B 


HcJ" 


5>1     j-1  o  tf  X1  *A     *  xi3  xi 

Let  us  indicate  these  new  coefficients  with  the  obvious  notation 


I(1)  =    a(l)  +  i^X 

Then  the  new  invariant  coefficients  (associated  with  D),  for  the  terms 

a* 

and  in  terms  of  the  original  coefficients,  these  "become 


>3ll  4-  .^«L  4f*(t)l-  *U»V)(  tj)  f    >1,, 


*XJ 


9? 

To  show  that  3110  invariant  coefficient  for  terms  of  lesser  order  is  a  true 
Invariant,  we  must  first  investigate  further  the  character  of  these  other 


S<  We  will  now  give  a  general  method  for  the  determination  of  the  in-> 
variant  coefficient  of  any  order  term  in  an  arbitrary  remainder  Be 
Consider  any  identity  given  by  a  sequence  of  operators* 

D(a)    »W      <  ...  D(q)  vDiS)  ,«*»«  +  *■ 

V*#1a  V*  fe    V   q.  ±#   « 

where     (n-a)  -I-  (n-b)+-  •••  +  £n-s)  as  n. 
We  will  denote  for  "brevity 

D(a)  d(<i)      »£)r 

Suppose  we  wish  to  find  the  invariant  coefficient  for  the  derivative 
u^'  .  when  e  »       is  a  term  in  the  expression  J>\         1   (u). 

S$nee  each  differentiation  ^  ,  i  =  l,  «».n,  never  appears  in  more  than 
one  operator  on  the  left  hand  side  of  (15),  we  see  that  the  left  hand 
side  of  (15)  contains  the  term 

£^0)  *(P)    n,   , 

and  no  other  term  on  the  left  hand  side  of  (15)  contains  v^'  as 

a  factor.  Thus,  recalling  the  procedure  necessary  to  make  D(u)  appear 
on  the  right  hand  side  of  (15)>  we  see  that  the  Invariant  coefficient 


for  u       must  "be 

Mil  .o.ffi. 


(16)       *C(e0)  -  ami< 
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Now  suppose  we  wish  to  find  the  invariant  coefficient  for  the 

(p)  ( s}   /  \ 

derivative  u  '    >  when  there  is  no  term  in  the  expression  3r  . '  .  (u) 

which  contains  u^p'      as  a  factor.  We  may  assume  that  D,  '   i  (u)_ 
m^c.m  ^...jg 

(p-j) 
contains  a  term  fl0  ^   0m     for  seme  J=l,2,  ...p,  where  we  set 

(o) 
u*  '=  u.  This  statement  is  ©Tbviously  true  for,  at  worst,  the  expression 

^o 

(  a) 

V*  •   (u)     must  contain  a  term    eQu,    (although  e     x&Y  vanish  identically). 

1     *ls 

We  may  then  apply  LeRoux's  rule,   extended  to  an  operator  with    n    inde- 
pendent variables,  to  the  product     e J*  •  -   Since  the  original 

*    °  m. . .  .m.     . 

J.  P"^  X 

operator    D    never  contains  any  differentiation     a     multiply,  for  any 

1  •1,2,...,  n,  the  same  is  true  of  the  operator  JJ  „     Thus  our  extension  of 
LeRoux's  rule  for  this  type  of  operator  is 

(175     D(uv)  =  uD(v)+T       |u    d'    (*)+l      JT     J£*        p*  (v)      .„., 

and  hence 

&-^-«3) 

£=1  p-J+k 

~r  ...  t  os-7^  „  (c  )u    .» -_  *  •  ••     • 

■k-d+itr:1^  °  V"mp 

(p) 

From  this  we  easily  conclude  that  the  invariant  coefficient  for    ul 

must  he 

(J) 
(19)  ^  (c0)   -  a 

V,HlooomP  V0°mp 


9k 

We  must  point  out  that  formula  (19)  does  not  give  the  complete 

picture.     For  If  Jj  is  the  identity  operator,  then  it  must  "be 

VWomP  /_% 

that     c0  — a^  ^^       ,  and  that     a  Bj*|     „       appears  on  the  left  hand  side 

of  (15) .     That  Is  to  say,   there  is  no  need  to  add  and  subtract 


a™        ™  u™        -     *o  the  right  hand  side  of  (15)*  hence  there  is  no  term 
ml0..mp  mlt- 


•  o  «.l 


(p) 

containing  %  ,ooII.  in  B.  In  other  words,  there  will  be  no  invariant 

coefficient  for  vS^'         •   The  same  conclusion  holds  whenever 

m-j_...ia 

(a) 

a       appears  as  a  coefficient  in  any  one  of  the  operators  D     .  , 
m-|_o.»mp  *l*#*  a 

D;j  '    •  ••>  D/.     •  This  is  apparent  from  the  left  hand  side  equation 

(11),  and  the  discussion  which  follows  it  regarding  the  determination  of  each 

term  in  D(u)  which  appears  on  the  right  hand  side  of  (ll)«.  To  summarize, 

whenever  a    -.      appears  as  a  coefficient  in  any  operator  D»a'   .  , 
^l***  p  l*'*  a 

(h)  (s)  ,   v 

T):   *      .     ,  ...  Di    i  ,  which  is  on  the  left  hand  side  of  (15),  there 

will  he  no  invariant  coefficient  for  u^'     appearing  in  B  on  the 

V'^p 
right  hand  side  of  (15).  Otherwise,  there  will  he  an  invariant  coefficient 


may  he  zero)  appearing  in  B  on  the  right  hand  side  of  (15)>  which 
will  he  given  hy  formula  (16)  or  formula  (19)  as  appropriate. 

F.  As  yet  we  have  said  nothing  about  the  character  or  number  of  quasi- 
invariant  s,  i.e.  the  invariant  coefficients  for  derivatives  of  u  of  order 
less  than  n-£.  As  we -have  pointed  out,  however,  there  will  he  no  invariant 

coefficient  for  a  derivative  u^p'     when  the, coefficient  a       appears 

m*  •  o  0m.  m-^ . .  .m_ 

on  the  left  hand  side  of  any  such  identity  (15)«  From  this  fact  we  may 
easily  compute  the  maximum  number  of  quasi -Invariants  associated  with  any 
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Identity  (15). 


total  number  ©f  terms  in  D(u)  which  are  of  order  less  than 

n 
n-2,  provided  no  coefficient  is  zero,  is  simply  ^T  fe\  .  If  we  sub- 

tract  from  this  the  total  number  of  coefficients  of  derivatives  of  order 
less  than  -n-2  in  D(u),  which  appear  on  the  left  hand  side  of  (15), 
assuming  all  are  not  zero,  we  will-have  the  maximum  number  of  quasi- 
Invariant  6/ which  can  appear  in  R,  on  the  right  hand  side  of  (15)  •  For 
if  such  a  coefficient  would  appear  on  the  left  hand  side  of  (15),  were 
it  not  zero,  then,  regardless  of  the  fact  that  it  is  zero,  there  will  he 
no  quasi -invariant  corresponding  to  this  term.  This  number,  then,  is 


&        "  n-a       n-b    x      n-s 
1*3  '  1-3      1=3    '      1*3 


(n-a)f  (n-b)*-  ...  f  (n-s)  -  n,  where  if  c<3,  we  set  ,^Jftrr( 


10. 
For  n  2,  the  number  of  quasi -invariants  is  obviously  equal  to 


zero.  For  n»3>  we  observe  that  an  unique  quasi -invariant  appears  as  the 
coefficient  of  u  in  each  identity,  hence  there  are  a  total  of  12  quasi - 
invariants  for  the  entire  third  order  case. 

For  n>3,  however,  the  picture  becomes  somewhat  clouded.  This 
occurs  because  while  no  quasi -invariant  can  appear  twice  in  the  same 
identity,  the  same  quasi -invariant  can  and  does  appear  In  more  than  one 
identity.  Thus  the  determination  of  the  total  number  of  quasi -Invariants 
for  a  given  order  n  seems  to  be  a  rather  tedious,  lengthy,  and  unre- 
warding problem.  To  obtain  an  indication  of  the  enormity  of  the  task  in- 
volved, however,  it  is  quite  simple  to  obtain  upper  and  lower  bounds  on 
this  number,  which  we  shall  designate  M(n) „ 
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An  d^sHediate  lower  bound  Is  a  consequence  of  the  fact  that  am 
uniqu©  quasi  =lnvar last  results  in  ©ach  Identity  as  the  coefficient  of  u0 
Hence  N(n)  Is  a  lower  "bound  for  M(n)0  To  obtain  an  upper  bound,  we 
merely  observe  that  the  total  number  of  quasi -invariants  associated  with 
D  must  be  less  than  the  number  which  would  exist  if  ©very  quasi -invariant 
appeared  in  only  one  identity 0  Now  with  ©very  ordered  integral  partition 


of  y 


Ik 


{•an 


*£••***!  J  ^itf* 


there  are  associated 


(■)(■ 


^ll  °<  2  A   oc, 


r 

oc, 


identities  of  the  form. 


(21)     D^  B^^iy^LA^^u)  ^  D(u)  4-  B,   ( 

equations (6)  and  (7  ))o 

With  each  of  these  Identities  are  associated  a  maximum 
invariants,  given  by  (20),  namely 


see 


k   &<1 


(22) 


n       fn-J   ,  .    k   051 


Thus  if  we  sum  over  @\   ,  the  set  of  all  ordered  integral  partitions 

i 

of  J,  and  then  sum  over  J,  we  will  obtain  the  number  of  quasi^invariants 
which  would  exist  If  ©very  quasi-invariant  appeared  in  only  one  identity,, 
This  number  is 


(23) 


M  (n) 


n«l 


(;) 
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Then  for  n  £  3,  wo  have 


N(n)  ^  M(n)  ^  M*(n). 


#- 


Computing  a  few  values  for  M  (n)  w»  find  the  following: 


n 

N(n) 

M  (n) 

3 

12 

12 

k 

7^ 

362 

5 

54o 

'8510    . 

G.  Now  we  are  in  a  position  to  complete  the  proof  of  Theorem  V,  We  must 

show  that  no  invariant  coefficient  for  terms  of  order  less  than  n-2  can 

in  general  he  a  true  invariant.  We  shall  prove  this  rigorously  only  for 

identities  in  which /j—  9__  ■+■  a(  \t     an<3-  shall  only  indicate  the  proof 

for  other  identities. 

Consider  the  invariant  coefficient  for  u       ,  where  p<*-2, 

p  m^.. nip 

and    m .    afc  r     for  any     i  —  l,2,oe.,pe     Then     2-  c.u  is  contained 

1  iaQ     1  m1...mp_i 

in    D  (u)     and  in  fact 
r 


(25) 


co 

a. 

Vi 

cl 

- 

a 
m. ..« 

Vl*> 

CP 

=■ 

a 
r   ) 

where  of  course  r  must  "be  properly  ordered  with  m,,  ...jEt. 
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Now  consider  the  change  ©f  variables    u  *  A    u,  where  \-\(*)  f=   Oa 
We  extend  the  equations  (]A).to  obtain 


all  *\  J1«P       v 


o        u 


(26)     &  xx°"*  ^Xi°l^  Xi4I"  ^Xn  ^xm  •••^ 


n«»p«l 

>       X       ,     1 

pU    p«n=l        P+i^l        a 


•        »       • 


o  #  i  =,  n»p 


aptt 


X 


s    AU! 


nip 


V^p 


+  .. 


If  we  substitute  (26)  into    D(u)  »    0,  and  divide  through  by    A    , 
we  obtain    D^l).     The  coefficients  in  which  we  are  interested  are  then 

n 


(27) 


n 


4^>  a  1      O  X         +■ 

j,k»p+.l      x        s  J  K  A    rfm1 

v  n-p  . 
4...f    1      O        A 


■»  *>  t»  *  *         00 


?      > 
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a™ m  r  =0^   _  -  r      ■     '"£""     a*,         ■»  ™  r       1    o  A 


J.pfl  A        mj 


J*G 


n 


j,icp*i   ^  -vr*    A  ax  ^  T 

j*k 

1      ^  A 

■4— r    -^ L2 — -. — _.      ,     where     z      =.x       ; 

and    ¥(rj   *    a^  ^.  _1    ^X    . 

In  all  the  above  we  nmst  comment  that  the  subscripts  involved  must  all 

he  properly  ordered  in  accordance  with  the  conditions  of  (5). 

The  quasi -invariant  for     u  m      in  this  particular  identity  is 

!••     p 


simply    Q  =  /T(c0)   -    ^.j.    ,  or 


a 


Q=  a     (a  )  +  a,   %(a  )   -  * 

^— —      m,,»»m  r  \v)     a,..jnr  m-.e.m-  m 

Then  the  transformed  quasi -invariant,   for     u  ,  will  he 

m-i  •  •  •mY> 

(28)        "5  =      ^     (a  m  J  •*   -       (a  )   -    T 

v      '         ^        -c-~—      m  • .  .m  r'       a/    n v  hl.  <>  •  oia  r'  m_  •  •  .m 

3x~      ip  \*)    '"i      p  "i      p  • 

Substituting  relations  (27)  9   we  find,  after  a  somewhat  tedious  computation 
that 
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n 


J#p+^ 


A  7T 


J,k-P4l      r  *  ^j*  ■* 

j*k 


(29) 


+    .    .    .  + 


n 


J-P4i 
J*P4^ 


&< 


a(4)'  +  a(r)a(j) 


-  a         )  O  A 

mp+l        J-l   mJ4>l-      "P-**     "»4<a. 


We  observe  that  the  coefficient  of  eacli  derivative  of  A  In  this 
expression  is  itself  a  quasi -invar iant  which  arises  from  the  particular 
identity  we  have  chosen. 

Then  using  an  obvious  notation  we  write     (29)     as 


'  *  Q    +  X  Qj  -r  -f. 


X  5^ 


(30) 


^  (n-(p+2)) 


(jr)  ■  (jr) 


Now  if    p-^n-2^ 


Q   *  Q   +■  Q~  1  i\ 
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But  Q  sO   (a.    )  +■  a/  \a„     -a.,   /-%/,%    — 
t»  <rr   l...n'  '  a{r)   l.o.n    l.<>.(r-l)(r+l)  .«.n 

^  U)  +  *(r)  -1  "  a(r)  *  °>  W° 
w  r 

Q  ■  Q  as  was  proved  previously  in  the  more  general  case. 
If  p>n-2,  however,  then  Qr  Q  -fa  linear  combination  of  other  invar- 
iants arising  from  the  same  identity.  Thus  ^^Q  in  general,  for  the 

particular  identities  D„    _   _    *  D  (u)  —  D(u)  +  R. 

J-   r-l  ttl  ft  ^r 

Although  we  shall  not  prove  it  rigorously  here,  it  easy  to  see  that 
for  p>n-2,  no  matter  which  identity  we  choose,  therms  involving  o  A   , 

^ ,   LJ ,  . ..,  will  arise  in  the  expression  for  any  transformed  quasi- 

invariant,  since  these  factors  appear  in  the  transformed  equation  D(u)=0, 
Thus,  we  assert  that  for  any  quasi -invariant  Q,  Q*  j*  Q,  and  in  fact  we  may 
expect  that  QmQ,  +  a  linear  combination  of  other  invariants  arising  from 
the  same  identity.  This  is  the  justification  for  the  name  "quasi-invariant, 
and  this  completes  the  proof  of  Theorem  V. 

H.  In  order  to  see  the  problems  involved  in  applying  the  cascade  method 
for  the  n  k  order  operators  we  have  "been  considering,  let  us  consider  a 
particular  kind  of  identity  once  again.  We  will  study  identities  of  the 
form 

(3D    \  4-U-i  wi..*w  ■  \  (^fa(r)u))= 
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n  \  ^— 

-D(u)  4-  J_    /d^£j+.(p)*(4)-«{jr)  \  *{3r)  '♦iSj.H^VlH 

[3D  *r         3  *J 

viler©     Q,  .  is  the  quasi-invariant  associated  with    u.  .         , 

ipfr,     p=l,...n-3. 

Denote    d}21"1'  ,     n         (u)  =  3  u  +a/   \u  =  u*       T    « 
l...r-lr+l...n        -\„     ~  W  >       J 

—  ^a(r)  "H^   ^  a/a   "  aM  ^  '>■■     aa^-  observing  "tiia-t    D(u)=sO,  we  write 
(32)  S2£         <r>  •> 


V 


'  (u')  =    T"     T-U/  ^\"i"2lQi  1        U4  i 

j3         (J  T.--#1a-3  11,'#  n-3* 

It  would  "be  most  desirable  if  all  the  T's  and  all  the  Q's  in  (32)  were 
identically  zero,  for  if  this  were  true  we  would  hare  succeeded  in  reducing 
the  order  of  the  original  equation  "by  one,  resulting  in  a  single  first 
order  equation,  coupled  with  an  equation  of  order  (n-l).  Ideally  "by  this 
method  we  would  like  to  reduce  the  nth  order  equation  to  n  first  order 
equations,  hut  this  happy  result  will,  unfortunately,  occur  very  seldom. 
In  the  more  probable  event  that  not  all  the  T's  and  Q*s  vanish 
identically,  we  would  then  hope  at  least  to  he  ahle  to  cascade  the  equa- 
tions in  the  manner  of  the  preceding  chapter py  until  such  time  as  the  T*s 
and  the  0,'s  may  all  he  identically  zero.  To  this  end  it  will  obviously 
he  necessary  to  impose  certain  rather  rigid  conditions  upon  the  coefficients 
of  D(u)  •  Let  the  hypotheses  he 
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(a)  T,=  ...T  ST  .S  „,  =  T  rQ.   ,.M    =T#  0, 

1     M  r*l        n   i1    i^     t- 

(b)  all  coefficients  are  f mictions  of  X-  only, 

(c)  all  coefficients  with  an  "rM  subscript  are  identical. 


Then 
(33) 


Dr(u')aT 


u-f    i 
1    n-2 


,  if*. 


Solving  u  =  d  u  +  a,  ,u  ,  for  u.  we  find 


u  =  e 
hence 
u 


K)<*r   £  j.for)**  ,.  ^^...^  ^  ...^ 


VV.2 


.3'      t, 


-K 


Ndx 
rj  r 


0  *<  o  •  •  £>■*-■» 


J 


far) 


dx 


u  dx^. 


f  F 


•"n-2 


* 


a,  xdx 
(r)  r 


-  e 


xi  ...<Jxi  (J         r   J 

xl  n-2 


where  q  is  the  nuaber  of  non-zero  i  • 

P 

Substituting  these  expressions  into     (33)     we  obtain 


for)*** 


»       .       I 


»P  (»  )  «■ 


Me  ud^  +  Fj, 


We  then  take  the  partial  derivative  of  both  sides  of  this  expression  with 
respect  to    x„,     to  obtain 


Ca,    vdx 
0   (*0     r 


a(r)Dr(u*)  +  J^UCU1)    .     ^logs  iC(u')   = 


j^Cr)^ 


5<1     . 

u 


H+r  . 


lOii 


Cancelling  out  the  exponentials  and  collecting  terms  we  find  that  we  have 
an  expression 

jD  (u»)  -    0  , 

th 

where  ,D  is  an  n       order  linear  operator  of  the  form  (5) $  but  with 
coefficients  differing  from    D.     These   coefficients  we  may  compute  to  be 

e'/tO    =:    a/_^-  "^logT        _ 


(r)    -    a(r)' 


2v 


a8      ...     —  a-    <,.„„  where     i     =r     for  some     p-l^-««n5 

il       in  1       xn  p 

*\...l        -^ai!...i      r  %       i     /•*•         -T,  where 

±         n-2  ~    J         n-2       +         1...  n~2  (r) 


i  ^r  for  any    p=l,  ...5  n-2   0 
P 

If  now  we  compute  the  new  invariants  we  find 


T  •  =  T  -     jLIli 


and     Q» „  „  ■  y     - 

xl  •  •  oln-3 


rv 


But  by  hypothesis     (c)}  &( -,)     -  a.  r    for  all  j^.  r    and  all  permuta- 

tions  ^•••in.'3   5  with  i  •£  r     for  any    p=l>    •  «.*  n-3» 

Thus  if  the  new  invariants  are  all  zero 9  we  may  reduce  the  order  of  the 

equation    D(u')   =  0    by  one.     If  the  new  invariants  are  not  all  zero5  the 

hypotheses  of  this  method  are  again  satisfied.,   and  we  may  cascade    the 

equations   in  the  hope  that  the  invariants  will  eventually  become  all  zero. 

This   can  only  happen  if     T  =  n     ^J^jA)   s  where       n       is  some  positive  integer, 

3  xr 
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Other  methods  can  he  developed  for  Identities  of  this  form  and  for 
identities  of  other  forms,  hut  in  every  case,  severe  restrictions  must  "be 
placed  on  the  coefficients  hef ore  cascading  may  h*gin.  We  see  then  that 
the  prohlem  of  attacking  these  equations  in  this  manner  "becomes  extremely 
difficult  and  restrictive  as  n  increases.  For  n>k     the  prohlem  "becomes 
enormous,  with  staggering  numbers  of  invariants  and  identities  to  he  manipul- 
ated, and  a  large  number  of  restrictions  on  the  coefficients  necessary  to 
cascade  the  equations.  If  it  is  desired  to  cascade  the  equations  it  is 
recommended  that  only  first  order^  substitutions  he  employed. 


SECTION  71 

\ 

A  BRIEF  SUMMARY  OF  EARLIER  EXTENSIONS 


A.  Extension  of  the  Laplace  cascade  method  have  "been  made  by  many  not- 
able mathematicians.  One  of  the  earliest  of  such  extensions  was  made 

21 
"by  Darboux  himself.  '   Darboux  considers  the  following  system  of  equa- 
tions of  second  order,  with  n  independent  variables. 

Let  poj  Y\j  Pp.-*  '  '  'J  £v»  k®  a  s3rs'i'®m  °^  n  independent  variables, 
and  consider  the  system  of  equations    , 

(1)  %  |ia  =  aUD^.+a:klii  '»   ^l»0»l#2,...n-l;  i*k  # 

This  system  contains  n(n-l)  linear  hyperbolic  equations  of  second  order, 

2 
and  we  observe  that  this  is  an  overdetermined  system,  whereas  the  system 

considered  by  us  in  Section  III  is  exactly  determined.  Darboux  seeks 

first  a  necessary  condition  to  give  n  linearly  independent  solutions 

in  addition  to  the  trivial  solution  u  m  constant.  To  accomplish  this, 

Darboux  forms  the  third  derivative  by  taking  &■   of  both  sides  of  (l). 

Interchanging  the  indices  1  and  k,  and  equating  coefficients  of  like 
derivatives  of  u,  we  obtain  the  relations 

(2)  \i&   -   ai  a  k  +  aliaik  "  aik*0Jc  (l«fe*i). 

Interchanging  the  indices  i  and  1  does  not  change  the  right  hand  side 
of  (2)  and  hence 

2aik  _  £alk 


^     ^i 
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Holding  k  constant,  we  e©e  from  integrability  conditions  that  there 
must  exist  a  function,  call  it  logE,  ,  such  that 


(3)  **~-±-*Jk.  (i*k) 

Hk  a^i  3 


Then  condition     (2)     "becomes 
and  the  system     (l)  has  the  form 

(5)         £»  -  i  ^  la  .  JL  3  Hi  .5  u ...  o 

l,  &  ■  u, x, <-,  •  o  o,n**i|  iriCo 
Using  condition  (k)     it  is  easy  to  verify  that 


(6)  fik  (H.)s  0  9       i#k*l. 

To  start  the  Laplace  method  we  assume  that  we  have  an  equation 
of  the  form  (5),  satisfying  the  conditions  of  integral) ility  (k)  „  Thus 
we  consider  the  particular  equation  (for  fixed  i,  k), 

(7)  fik(u)  -  0  # 

Darboux  then  considers  the  substitution  function  V  defined  by  the 
relation 
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(8)  Su  -  _l_^Hk(u  +i/). 

If  v©  substitute  (8)  into  (7),  the  equation  then  becomes 
l^k  3U  ,/dlo^T,  .  i  Jl   ,  ^  H.\  , 


We  see  that  our  invariant  h.„,   is  now 

ik 

and  if    Iw—  0,   the  equation     (9)     has  the  simple  solution 

V  s  constant  # 
We  may  gain  a  little  insight  into  these  proceedings  if  we  adopt 
the  notation  of    Section  Y.     Thus,   for  equation     (7), 

¥iJfk        Hk    2f i  ^k      Hi    <^k  3pL 

-    <2  -  a       <)       -  a,        O 

and  "    5^        lk  3fk      ^    <5pT  ' 

Then  the  substitution  function    V    is  defined  by 


(io)  Dikpk(u)  -  V*0-- 


vhile  tho  imrariaat    h..      ™W  t>«  -written 
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(n)  K„  ■  *>  '    CO, 

I 


k?J 


and  hence  (9)  "becomes 

(12)        %(D*Pk  (u))=  ^(u)+(D^1(-a^))u  • 

If   D    (-a  )=0,  we  see  that  (10)  reduces  the  order  of 
ikpv   ik     ' 

(12)  by  one,  and  the  equation  can  he  solved  by  quadratures  „ 

In  the  more  likely  event  that  D*       fas-J)   t^  0,  we  may  cascade  the 

ikpi  ik  I 

equations  in  the  following  manner  0  Let 


and 


Then 


5^    -    J*  (u^) 


and  u»  satisfies  the  equation 

(13)       ^V    -  JL  <^£  3V   -  JL.  ^Lj     'bV    -     0 

aft^     I*  aPi  apR     i-i    ^K  a^ 

In  like  manner,  if  we  introduce  the  quantities  L»  ,  defined  by  the 
relations 
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h,   £3L       ak|^  (k'*i  ^k) 

<^i  0?i  J 


then  the  function     ^  ,   far  distinct  values  of     i*     and    k',   satisfies 
the  system  of  equations 

^V3fV      V   ^fk'  *&  V  dpi*   ^k'  ^ 


is  of  the  same  form  as  (5).  •>  Hence  we  may  iterate  the  substituion, 
and  cascade  the  equations  until  such  time  as  the  corresponding  invariants 
may  vanish. 

(22) 

B.  In  l899>  Jo  Le  Rows.  extended  the  method  of  Laplace  to  linear 

partial  differential  equations  of  greater  order  than  second.  Le  Roux, 
however  considered  the  equation  with  only  two  independent  variables, 
while  we,  in  Section  V,  have  considered  the  equation  in  n  independent 
variables.  In  the  notation  of  Le  Roux,  the  equation  to  he  considered  is 

the  n^1  order  equation 


(15)     D(z)  =   Z-       ft!  A 

where  a^+£5#  n,  0<  f*  n,  (3^  n. 

If  we  suppose  that  the  highest  order  of  differentiation  with  respect  to 
x  in  equation  (15)  is  equal  to  n-  p,  then  the  collection  of  terms 

which  contain  such  a  differentiation  are 


+  ...  +  6  * 


n-p 


v     z 

$7^    ' 


Ill 


W©  denote  th©  expression 


o    +  h  d  "  j. 


ri-l 


o.  4-g 


as  the  differential  multiplier  of  the  term  O         z   ,     and  we  further 
assume  that  the  coefficient  a  is  equal  to  one.  L©  Roux  considers  the 
equation  (15)  from  the  point  of  view  that  there  exists  a  particular 
integral  of  the  form  of  Euler 


(16) 


(m)  .     (m-1) 
u  X   +  UiX 


.  +u 


& 


where  X  is  an  "arbitrary1  function  of  x,  and  the  coefficients  u 

are  functions  of  x  and  j» 

If  we  regard  D(z)  as  a  polynomial  in  &   and  j&  ,   we  nay  introduce 

5*     Jy 
the  notation  of  Section  Y.  We  designate  /_.  n\         as  tht  expression  ot>- 

XV 

tained  when  we  differentiate  D  with  respect  to  <7  9   p  times,  with 

Jx 
respect  to  ^d  3   q  times,  and  apply  this  operator  to  the  variable  Z.  For 

th 
instance^,  the  (p+q)    derivative  of  the  term 


is 


oc(*-l) 


3x*3y* 


d**P}y*-* 


and  hence 


(IT)  D     (z)^n(n-l),80(n-p. 


<-P<-(3-<! 


(n-p-q^        A  y*  ^  *  Z, 
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Since     (15)     is  assumed  to  have  order     (n°p)     with  respect  to    x,     we  see 
that  the  coefficients  of     (l6)     oust  satisfy  the  equations 


(n-p) 


a~TT  ^   (UJ   "   ° 


(18) 


r   1    _.      Djn-P       (v.)   +  1 

[n-P)  S  (n-p- 


(n-p-l) 


•       o       • 


In  the  case  when  the  variable  x  is  a  simple  characteristic  variable,  Le 
Roux  considers  the  transformation 

*i  =    i.  V1  (z)  s     f£-  +  ^21-1,0  z> 

and  defines  the  functions  ^  and    v/.    - 


(19) 


4" 


1  -z. 


jMn~l^   ,     U    -   ^ 


1- 


tt 


We  then  consider  the  equation  resulting  from  (15)  "by  means  of  the  trans- 
formation (19)5  and  collecting  on  the  right  all  terms  which  do  not  contain 
any  differentiation  with  respect  to  y,  we  obtain 

(20)  m^Y-K^  X^a.n,...^ 

where    /^  ( IA )     designates  a  differential  expression  of  order  (n-l),   in 

^n-1         11 -l 
which  the  coefficient  of      &      *i$x      "  is  equal  to  unity.    The  coefficients 

are  given  "by  the  formula, 

1  AnA     ,dy 

(21)  X     -     e0      n"X 


(P-D,      , 


Uq       IP' 


1      Dxp-i  (ti0)  , 
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The  order  p  of  the  right  hand  side  of  (20)  is,  in  general,  equal  to 
(n-2).  It  will  1)©  less  only  if 

(n-2) ,   x 

D  n-2  ( V  ~  °  • 

x 

Le  Roux  then  proceeds  to  define  the  invariants  h,  "by  the  follow - 
lag  rule.  If  the   X's  regarded  as  functions  of  y  are  all  linearly 
independent,  then 

JL  •  e  o         «L 


(22) 


1 


By 


0A1      o  ..dXj 


d  Ai  ...  g  Xi 

h7i     a71 


which  are  (pfl)  in  number,  setting  h  «X  ♦  If  the  A's,  regarded  as 
functions  of  y,  are  not  linearly  indepea&ant,  we  define  A^  as  the 
first    A    which  is  linearly  Independent  of     \Qt    Ap  as  the  f Ifibet  X --vhich 


is  linearly  independent  of  \^    and    A  -, ,   e 
is  replaced  by 


(22') 


h. 


T7 

0   •   • 

a  Ao 

a*1 


^AoCt        ^A 


:c.     Then  the  determinant     (22) 


0Ci 


I1 


a  •  a  o  •  o 

<}  A°a    ^  Xxi 


a/ 


a/ 


To  show  that  these  are  generalizations  of  the  Darboux  invariants,  we  must 
consider  the  transformation 

z  ■  z       f(x,y) 
(23)  |x"  <?(*') 


Ilk 


The  first  of  these  transformations  is  effected  without  changing  the  value 
of  V   ,   siace   u'  s     o    ,  and  hence 


z\  z 

o     o 

The  transformation  of  the  independent  variables,  however,  causes  equation 
(20)  to  "become 


%xi2  ^"1 


where 


(      ) 


if' 


,n-l-i 


3y'UV  2  X 


•      o      a 


the  coefficients  0  "being  dependent  only  on  x.  If  we  use  the  relations 
(25)  to  compute  the  new  invariants  h|  ,   we  will  find  that 

(n-2-p)(i+l) 


1  VJ7'  jtV  -3* 


Since  in  general  p»n-2,  the  ahove  relation  reduces  to 


h 


(1+1) (1*2) 

hi  =  /q?  2  x\ 


2     fc1. 


Thus  the  determinant  h„  is  reproduced  mutliplied  "by  the  (i+l)  (i+2) 
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power  of  the  Jacot>ian  of  the  transformation,,  This  justifies  the  name 
" invariant"  for  the  determinants  h.,  and  LeRoux  calls  these  the  first 
generalization  of  the  Darboux  invariants. 

Now  for  the  equation  (15)  to  admit  an  integral  of  the  form 
z  m,  UgX,  X  "being  an  arbitrary  function  of  x  only,  it  is  necessary 
and  sufficient  that  the  coefficients  A  of  (20)  are  all  zero0  This  is 
evident  if  we  but  consider 

D(uX)«  X(D(u)  +  xV(u  )4  x"D,V>(uJt...tX(n-l)  D^  (uo)~0  # 

1X21X  (n-1)  I   *" 

Since     X     is  arbitrary,   and  this  equation  must  be  identically  zero,   it 
follows  that 

which  proves  the  assertion.     Thus  we  see  that  the  vanishing  of  the  invar- 
iants plays  the  same  role  for     (15)     as  it  does  for  the  second  order 
equation. 

LeRoux  then  considers  the  case  when    x     is  a  multiple  character- 
istic variable,    so  that     (15)     has  the  form 

,    *      ^       ^  n~aL  ^       ^  n-a-1 

(26)       DU)  =  <p0    d z_+<Px    d    .       f +  ...  +  <P.«*0. 

N       n-m  *\ nna-i  n-a  •? 

d   X  (jX 

where  %j^fij  '  '  '  designate  the  differential  multipliers,  e.g. 

<$><,=    jjL  +  ai  |^r  +  aa   ^  +  ...      (r.*»)..  b 

The  first  substitutioa.  we  have  already  considered,  namely 


aT 


that 
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We  attempt  to  decompose     <¥0{z)   into  differential  factors,   ao 
<P0(*)»*l«2   •••  ^r4uol  i 


Then  we  define 


Jy^r-1  JF  «Cr^  "3y  tf."3jr 


y  <*j  ^y  uf 


n 


c 


y   u0   .) 


^lo  9    1    ^ 


y  k. 


_£  *  or  nU 
1  1  o 


3    1  *■  °SF^  -37  5^  #  1^  *i*a  \§j 


90t0«2«0  -> 


«         •  • 


amd  using  these  we  are  able  to  reduce  the  equation.     If    (P0(z)  does  act 
have  such  a  decomposition,  we  may  still  reduce  the  equations  "by  a  method 
which  will  be  described  "briefly.     The  proof  of  the  following  is  rather 
lengthy  and  detailed  and  does  not  hear  repeating  here.     It  may  he  found, 
however  in  LeRoux's    work  previously  noted  (footnote 22). 

Let     u  9  uit  ^V,,  •  •  •  denote     a  set  of  linearly  independent  integrals 
of  the  equations  (18)  which  define  the  conditions  on  the  coefficients  for 
the  existence  of  a  solution  of  the  form  of     (16) .     We  then  define  the 
functions 


kaJL 
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Set  i   _   v  k 

•  xi  m    :L  :i-2 

2 
ki°l 

The  equation 


h 


1  o      1   «oo  o     id 


J  li'^U-i      ?y  li  ^Uo)" 


admits  as  solutions  the  functions  u0u,.  «»ooUo.  Then  the  set  of  transfer* 
mations  is  easily  defined  "by  setting 


1        °  37 


(* 


s        u     1.       J      1    ^ 


Z2  "       ~o  J-l 


•  •  • 

z 


^     ^1,  L  —L^    JL  JL  JL-J: ^     id 


i-        ©*1A2 


°     JL*.     °         1     .„.    O      \    d    (z\ 


If  there  exists  a  particular  integral  of  the  form  of  Euler,  the  inrar- 
iants  1  will  become  zero  after  a  certain  number  of  iterations,  and 
the  chain  of  transformations  will  stop, 

23 
C.  Following  the  work  of  Darboax  and  LeRoux,  Dini   made  an  extension 

i 
of  the  Laplace  method  to  the  linear  second  order  equation  in  an  arbitrary 

number  of  independent  variables e  Dini  considered  the  equation 

2 


n         -  ~       » 


CT-i    iJ  3x73V    fin    J  J,: t    '  iJ      j1j 


(27) 

"1  .    jj  —  J.         ^  A  j 

and  posed  the  question,  can  we  transform  (27),  into  the  form 


n  ^<C}  n  -\ 

(28)     X  ki  ^~  +  M^  f  21  OC,  gl.   +L«  4-  H  =  0 

i=l  o   x±  i=i  i  ^x^          J 

where 
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n 


(29)  ©"     21    aJ  ^  41jz  J 

and  the  functions  a^,  ag,  ...,3^,  h,  k-^...!^,  K  are  (2n  +  2)  functions 
to  "be  chosen  arbitrarily?  If  w©  assume  that  this  has  "been  accomplished, 
substitute  (29)   into  (28),   and  equate  coefficients  with  the  left  side  of 

(27),  w©  olbtaim  the     (n+l)(n+2)    equations 

2 

(30)  V-j    4-    kjai-2Ait.       >  i>«]  =  l;-;*; 

n 

(31)  kjT*   +     a±M  4.^    kj^ai+oCi  -   Gi    j  i=l,...,n^ 

J.1    >x,  , 


(32)    Mb  +  y     k,  ^b_  f  I&I 
J-l    °   3 


for  the  (2n-f*2)  unknowns.  We  olb  serve,  however,  that  in  (30),  the  k„ 
and  a,  appear  only  as  products  kja  ,  and  not  all  of  them  are  zero  if 
(27)  is  to  "be  effectively  of  second  order.  Thus  we  may  select  arbitrar- 
ily one  k,  say  k^,,  and  set  L  ml,  without  affecting  these  relations j 
and  hence  there  are  actually  only  2n-l  unknowns  a^  and  k^.  If  we  consider 
first  the  equations  (3°)>  we  note  that  for  ns2,  these  are  3  equations 
in  3  unknowns,  and  this  system  is  determined.  For  n  £  3>  however,  the 

system  is  overdetermlned,  and  one  would  expect  to  have  (n-l)(n°2)  rel- 

2 

ations  among  the  coefficients  A^..  so  that  only  2n-l  of  the  equations  (30) 

will  he  independent. 

We  see  then,  that  hef  ore  we  may  even  "begin  to  seek  apply  a 

"Laplace  method"  to  equation  (27),  the  coefficients  Aj«  must  satisfy 

(n-l)(n-2)  conditions.  Dini  showed  that  in  general  these  conditions  may 

2 
"be  stated  as 
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(33)       V^s  "    *lA,'  «reel[v  -AhhArrVAhB  -  Wee  ,  V-g.3,,.^, 

j 

where  £  m  il,  £g  =  +1,  ami  these  signs  are  determined  "by  the  sign  of 

the  radicals  of  the  roots  of  the  quadratic  equation 

W"  SAlr^r+AIT=0- 
Here  A. .,   is  considered  to  he  the  first  noci-vanishing  coefficient 
Ajj,  i*l, .o.,n.  If  all  the  A.,  are  zero,  the  equations  (33)  are 
satisfied  Identically  for  all  Aj^,  ial;  0.o,n. 

The  classification  of  the  various  allowable  types  of  equation 
(27)  "becomes  the  next  important  consideration.  We  list  here  without 
proof  soma  of  the  more  important  properties  of  (27)  derived  "by  Dini. 
First,  if  n>2,  and  all  coefficients  of  (27)  are  real  valued,  then 
for  any  pair  of  variables  (xr,x0)  for  which  the  partial  differential 
equation  is  not  parabolic,  the  equation  must  always  be  of  the  same  type. 
That  is,  ignoring  all  pairs  of  variables  for  which  the  equation  is  para- 
bolic, and  considering  all  remaining  pairs  of  variables,  the  equation 
must  always  be  elliptic,  or  always  be  hyperbolic  type,  with  respect  to*«ach 
pair  of  variables.  Another  property,  true  even  for  complex -valued  coef- 
ficients.!, if  Aj^  f*  °>  and  the  partial  differential  equation  (27)  is 
parabolic  with  respect  to  the  pairs  (x^x  )  and  (x^,x  ),  then  it  will 
also  be  parabolic  with  respect  to  the  pair  (x^x  ).  These  properties 
follow  from  the  conditions  (33) • 

Dini  then  defines  the  following  terminology.  Equation  (27)  is 
said  to  be  of  parabolic  type  provided,  that  It  is  parbolle  with  respect  to 
every  pair  of  variables  (x^x^;  **  there  exists  at  least  one  pair  of 


120 


variables  (x  ,x  )  such  that  (27)  is  mot  parabolic  with  respect  to  [ 
(s^Xg)  then  (27)  is  said-  to  be  of  elliptic  (hyperbolic)  type  if  it  is 
elliptic  (hyperbolic)  with  respect  to  every  pair  of  variables  (x^Xg) 
for  which  it_is__giot  of  parabolic  type;  otherwise  (27)  is  said  to  "torn   of 
mixed  tj£Qo     Then  Dini  proves  that  if  (27)  is  not  of  parabolic  type,, 
then  either  there  exist  exactly  two  systems  of  values  of  k  and  a  which 
satisfy  (30),  or  there  exists  none,  while  if  (27)  is  of  parabolic,  type, 
there  exist  at  most  one  system  of  values  of  k  and  a  satisfying  (30)  • 
In  the  first  case,  on©  ©"btains  one,  system  of  values  from  the  other  merely 
"by  interchanging  the  k's  and  the  a  "a.  Denoting  these  values  by  (k,a) 
and  (a,k)  respectively,  we  may  call  them  conjugate  sets0 

Now  our  first  thought  in  attempting  to  reduce  (27)  to  the  inte- 
gration of  two  partial  differential  equations  of  first  order,  is  that 
perhaps  (28)  will  reduce  to  a  first  order  equation  in  ©  alone.  In  ad- 


dition to  the  iBr^-liBd^X    saaditioms  previously  noted  as  necessary  to 

2 

reduce  (27)  to  (28)  and  (29),  this  requires  n-1  further  conditions 

to  insure  that  erf".,****  <»..  =  of  s  L  =  0„  These  conditions  arise  from 
x   d  n 

the  (n«fl)  equations  (31)  and  (32)  which  determine  only  the  two  un= 
knowns  b  and  M„  Assuming  that  the  equation  is  not  of  parabolic  type 
and  hence,  without  loss  of  generality,  that  it  is  not  parabolic  with  res- 
pect to  the  pair  (x,x  ),  Dini  shows  that  the  n-1  conditions  can  be  formulated 
as 

(3*0  I  -  Mb  -  B—  0 
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(35) 


A 


11 


12 


Is 


21 


22 


2s 


where  A  = 
s 


n 


r-1 


Gl  "  Al 


2     2 


G 


—  0    ,  s  =-3>^>  o..,n 


ky,  Q   s  ,  and  B 


a 


r-1 


k   ^b 

r  =r— 
7)  x 


We  may  look  upon  L  and  the  C?C«  **  (2^)  as  Playing  the  role 
of  the  lM.yaria3a.ts.  Thus  to  say  ,that  the  conditions  (3*0  aad  (35) 
are  satisfied,  is  to  say  that  the  invariants  of  (27)  are  all  zero. 
If  we  denote  D  (u)  =  M^  4.  ^"  k  d  m  ,  these  invariants  aay  be  written 


as 


(36) 


J 


L  =  I  -  D  ("b) 


0^,  r= 


i  -  G  -  k^  -  D  (a±) 


In  order  to  start  the  chain  of  equations  when  these  invariants  are  not 
all  zero,  Dini  imposes  conditions  similar  to  those  imposed  "by  us  on  the 
third  order  equations  of  Section  IV.  For  equation  (27)  to  he  cascaded, 
Dini  requires  that  conditions  (35)  he  satisfied.  That  is  to  say,  the 

invariants  OC   are  &^-   zero>  "but  ^e  invariant  L  is  not.  Then 

i 

equation  (28)  may  he  solved-  for  z,  in  terms  of  ©  and  its  first 
partial  derivatives.  This  expression  is  substituted  into  (29)  and  a 
new  second  order  equation,  for  O  ,  is  obtained.  Not  only  will  this 
equation  he  of  the  same  form  as  (27),  but  th©  second  order  coefficients 


A. .  are  reproduced.  Hence  our  new  equation  is 


122 


(37) 


i,J»l 


S2e   + 


The  new  coefficients  are  given  "by 


G  s  G  +  K   -  A  -k  L 

s      s      s      e   s  a  J 


(38) 


^a"  La^ 


^k, 


where  K  =  >.  a  <*  s  ,  and  H  ,  La  are  given 


B 


a 


r=l 


Since  the  relations  (30)  used  to  determine  (k,a)  for  (27)  involve 
only  the  coefficients  A,  ,9   it  immediately  follows  that  if  there  exists 
a  set  (k^a)  which  splits  (27)  into  (28)  and  (29) $  the  same  set 
(k,a)  will  suffice  to  split  (37)  analogously.  It  will  he  necessary 
however  to  compute  a  new  M  and  h  f   and  a  new  set  of  invariants  L  9 
and  ^C.  .  The  conditions  under  which  the  ^T,  are  zero  "become 


(39) 
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which  are  quite  similar  to  the  conditions  (35).  If  these  conditions 
are  satisfied,  and  further  L  =  W  -Mb  -B  =0,  then  (37)  reduces 
to  two  first  order  equations  as  desired.  If  L'  is  not  zero,  however 
we  must  again  require  that  conditions  (39)  he  satisfied  "before  we  can 
iterate  the  process  to  obtain  a  third  second  order  equation  of  the  form 
(27)  o 

Dini's  very  remarkable  result,  however,  is  the  following,,  If  the 
conditions  (35)  and  the  conditions  (39)  are  all  satisfied,  then  we 
need  impose  no  further  condtions  to  iterate  the  process  indefinitely. 
That  is  to  say,  if  the  invariants  °C  ±     and  OC  ^  are  all  zero,  then  all 
invariants  ^C±        will  he  zero  for  any  positive  number  m  of  iterations. 

This  result  "becomes  evident  if  we  compute,  say,  G   "by  the  rule  given 

s 
«  ji        ti  »» 

in  (38)  for  G  ,  substitute  G1  -  A  ,  Gn  -  A    and  G  -  A   for  the 

corresponding  terms  in  (39)  >  and  employ  the  rule  for  evaluating  a  deter- 
minant when  one  of  the  columns  is  a  sum  of  two  or  more  columns.  Each  of 
the  resulting  determinants  in  the  sum  must  vanish  as  a  direct  consequence 
of  conditions  (35)  and  (39) . 

In  summary  then,  in  order  for  the  equation   (27)  to  "be  cascaded 
as  many  times  as  necessary  for  the  L  invariant  to  vanish,  there  must  be 
satisfied  a  total  of 


conditions  on  the  coefficients,  The  first  (n-l)(n-2)  conditions  are 

2 
necessary  to  reduce  (27)  to  two  first  order  equations,  while  the  remain- 
ing 2(n-2)  are  necessary  to  permit  the  cascading  to  continue.  We  observe 
that  for  n*=  2,  no  conditions  are  necessary,  and  this  is  the  original  case 


22k 


considered  by  Laplace,  Legendre  and  Darboux. 

(2k) 

D,     Filially,  Mention  should  "be  BEade  of  the  work  of  Burgatti  who 

considered  the  equation  of  elliptic  type, 

(kO)  d    *.    3z_4  a  3z  +  b  ^_z  +  cz=0. 

dx2       Ty2  Jx  £y 

Burgatti  found  that  the  two  expressions 


E    = 

3a 
J7 

- 

9b 
J*' 

K  s 

0  a 
£  x 

+ 

37  + 

a2  + 

b2 

2 

^  _  ^ _     *s,     2  .  ,2 

-  2c 

are  invariant  relative  to  the  transformation  z  »  X  z  •  Further,  if 
H  and  K  are  "both  zero,  equation  (kO)     may  "be  reduced  to 


2  •     n2  • 

+  3   z   -  0  , 


5 

which  is  Laplace's  equation ;  If  H  is  zero  hut  K  is  not  zero,  then 


{kO)     takes  the  form 

while  if  K  is  zero,  hut  H  is  not  zero,  {kO)     becomes 

<?  xc         3  y*      * 

where  (X  and  S  are  functions  of  x  and  y. 
In  the  notation  of  Section  V,  let 


D  -   g   4  L_  +  ajL  -I-  b  5   +  c, 
ax2    ^y2      3x      J7 
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then  J)£    ■     2   o        4-      a 


D '     =     2    ^        +     h 

y  57 


and  the  invariants     (kl)     then  "become 


(te) 


H     -         Dy(a)    -  l£  (D)     „    (Dy(a)    -  c)    -  (%(*)    -  *) 
-.  ^  _. 

K      =         D  (a)  +  d'OO    -   2C  =  (l/(a)    -  c)4  (D  (h)    -  c) 

2~      "  2  2 

Since  the  change  of  variable  and  the  vanishing  of  invariants  does 
not  reduce  the  equation  to  two  first  order  equations,  no  consideration 
can  "be  given  to  cascading  the  equations  of  elliptic  type  utilizing  the 
invariants  H  and  K.  The  work  of  Dini  shows  us,  however,  that  this 
equation  can  he  cascaded  if  we  admit  the  possibility  of  complex  coeffi* 
c lent So 
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